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Summary 

An asymptotic description is derived for the interaction between a 
shock wave and a turbulent boundary layer in transonic flow, for a particular 
limiting case. The dimensionless difference between the external- flow 
velocity and critical sound speed is taken to be much smaller than one, but 
large in comparison with the dimensionless friction velocity. The basic 
results are derived for a flat plate, and corrections for longitudinal wall 
curvature and for flow in a circular pipe are also shown. In Part I, solu- 
tions are given for the wall pressure distribution and the shape of the shock 
wave. In Part II, solutions for the wall shear stress are obtained, and a 
criterion for incipient separation is derived. Part III contains simplified 
solutions for both the wall pressure and skin-friction distributions in the 
interaction region; these results are presented in a form suitable for use 
in computer programs. 


V 



INTERACTION BETWEEN A NORMAL SHOCK WAVE 
AND A TURBULENT BOUNDARY LAYER 
AT HIGH TRANSONIC SPEEDS 

Part I - Pressure Distribution 


A. F. Messiter 




1. Introduction 


In several recent studies, asymptotic methods have been used success- 
fully for the derivation of rational approximations which describe the inter- 
action of a turbulent boundary layer and a weak, stationary, normal shock 
wave. It appears that correct limiting forms of the equations can be deter- 
mined, that numerical or analytical solutions to these equations are ob- 
tained easily enough to be of practical interest, and that numerical accu- 
racy may be adequate for important parameter ranges. In the limiting 
case to be considered here, still for an unseparated boundary layer, 
the shock wave extends close to the wall, the upstream influence is 
small, and analytical solutions can be obtained for most of the flow 
field. Pressure distributions are derived in Part I; the wall shear 
stress and the possibility of predicting separation will be discussed 
in Part II. 

In many transonic flows of interest, there occurs a shock wave which, 
in an invi sc id-flow approximation, is normal to a solid boundary, at values 
of the Reynolds number large enough that the boundary layer along the wall 
is fully turbulent. Since the strength of the shock wave must decrease to 
zero in the supersonic part of the boundary layer, there can be no discon- 
tinuity in the pressure at the wall. It is observed that the shock wave 
becomes slightly curved and is displaced slightly in the upstream direc- 
tion. As the Mach number upstream is increased, still below the value 
required for separation, the shock wave extends further into the boun- 
dary layer; experimental results [l, 2] show an initially rapid rise 
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in the wall pressure, followed by a gradual decrease in the pressure 
gradient over a distance several times larger than the boundary layer 
thickness. 

Asymptotic descriptions of these flows, in the limit of infinite Reynolds 
number, have been discussed in References [3] through [lO]; in particular. 
Ref. [8] contains the first steps of the present work. In each of these 
studies, the representation of the undisturbed boundary layer in terms of a 
velocity-defect layer and a wall layer [ll, 12, 13] is regarded as providing an 
asymptotic description as the Reynolds number tends to infinity [14-18]. 

The pressure gradient in the boundary layer is large near the shock wave, 
and consequently the forces resulting from changes in the Reynolds stresses 
are of higher order than terms retained, in most of the boundary layer. 

Thus, as for laminar flow [l9, 20,21] , an asymptotic description of the 
changes in the mean flow can be obtained with the use of inviscid-flow equa- 
tions for most of the boundary layer. 

The form of the velocity profile, however, implies two important 
differences from the laminar-flow case. First, for an unseparated turbu- 
lent boundary layer the wall layer is extremely thin, and the displacement 
effect resulting from deceleration of fluid close to the wall remains very 
small, even in a large pressure gradient. Thus, if the undisturbed velocity 
profile is known outside the wall layer, an approximation to the pressure 
can be found without knowledge of the flow details near the wall and there- 
fore without any further assumption about the nature of the turbulent 
stresses. Second, for a slightly supersonic external flow the sonic line is 
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located at ctn arbitrary position (outside the wall layer) in the undisturbed 
boundary layer, depending on the relative sizes of the nondimensional fric- 
tion velocity cind the nondimensional difference between the fluid velocity 
and the critical sound speed in the external flow. As the Reynolds number 
tends to infinity, one can then study three cases, such that the ratio of 
these parameters tends to infinity, remains constant, or approaches zero. 

Adamson and Feo [3] considered an incident oblique shock wave in a 
flow with velocity only slightly greater than the sound speed, such that the 
sonic line is located very close to the edge of the boundary layer. The 
corresponding asymptotic formulation was shown to lead to a local- inter- 
action problem requiring solution of the transonic small-disturbance equa- 
tions for the local perturbations in the external flow, expressed in appro- 
priately scaled variables. The influence of the boundary layer is repre- 
sented on this scale through an effective wall boundary condition specifying 
a linear relationship between the streamline slope and the pressure gradi- 
ent. Melnik and Grossman [4] studied a normal shock wave having strength, 
as measured by the nondimensional pressure jump, of the same order as 
the friction velocity, so that in the limit the sonic line is at an arbitrary 
location in the boundary layer. Numerical solutions of the transonic small- 
disturbance equations were obtained for perturbations in the defect portion 
of the boundary layer and in the neighboring external flow. Changes in the 
wall layer were also discussed in each of these papers. Melnik and Grossman 
later [5, 6] obtained additional numerical solutions for axisymmetric pipe 
flow. At higher upstream speeds, which might be characterized as "high 
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transonic speeds," the shock wave is stronger but ^'he boundary layer can 
remain unseparated. For this case, a first approximation for the flow per- 
turbations outside the wall layer was given by Adamson and Messiter [8] , 

The shock-wave strength, although still small, was taken to be large in com- 
parison with the nondimensional friction velocity, so that in the undisturbed 
boundary layer the distance from the sonic line to the wall is much smaller 
than the boundary-layer thickness. The corresponding problem has also 
been discussed for an incident oblique shock wave [7, 9] . A brief prelimi- 
nary description of some of the present results was given in Ref. [lO] ; a 
few details have since been modified. 

In the present work, analytical solutions are derived which incorporate 
additional physical effects as higher-order terms for the cas^ first dis- 
cussed rather briefly in Ref. [8] , when the sonic line is very close to the 
wall. The functional form used for the undisturbed velocity profile is des- 
cribed in Section 2, to indicate how various parameters will be calculated 
for later comparison with experiment. The basic solutions for the pressure 
distribution are derived in Section 3. In Section 4 corrections are added 
for flow along a wall having longitudinal curvature and for flow in a circular 
pipe, and comparisons with available experimental data are shown. The 
restriction to weak shock waves is removed in Appendix A, and it is verified 
there that the simpler solutions of Section 3 are adequate. Some additional 
results pertinent to the asymptotic matching of solutions in the region of most 
rapid pressure rise, near the beginning of the interaction, are derived in 
Appendix B. 
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2, Undistarbed Velocity Profile 


Nondimensional rectcingalar coordinates X and Y are measured along 
and normal to the wall, respectively, with Y = 0 at the wall and X = 0 at 
some point on the shock wave, e. g. , at the intersection of the shock wave 
with the edge of the boundary layer as defined below. The reference length 
is a geometric length such as the length of the boundary layer from a lead- 
ing edge up to the shock wave. The nondimensional mean-velocity compo- 
nents U and V, referred te the critical sound speed in the external flow, 
are in the X and Y directions respectively, and the term p' V‘/p has been 
included in V. Here primes denote fluctuations about the mean, and p ' V 
denotes an average value. The nondimensional mean pressure P, density 
p, temperature T, and viscosity coefficient p are referred to the criti- 
cal values of pressure, density, and temperature, and the corresponding 
viscosity coefficient, in the flow just outside the boundary layer and ahead of 
the shock wave. The sum of the nondimensional Reynolds stress Etnd viscous 
stress, in the boundary-layer approximation, is denoted by t, and has been 
made nondimensional with twice the dynamic pressure, in terms of the same 
reference quantities. For later convenience the friction velocity u_^ is made 
nondimensional using the external-flow density: 

2 ”*"w 1 2 

u = = V U c. , U = 1 + c (2.1) 

Tp 2ef' e ' 

e 

where the subscripts e and w indicate values in the external flow and at the 
wall, respectively, and c^ is the undisturbed value of the skin friction 
coefficient, referred as usual to the dynamic pressure in the external flow. 
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The nondimensional difference between the fluid velocity and the critical 
sound speed in the external flow is € , and in the present case u « e « 1. 

T 

For simplicity, an adiabatic wall is assumed and the total enthalpy is taken 
to be uniform. The ratio of specific heats is y and is constant. 


As in references cited above, it is assumed that the undisturbed 
boundary layer can be described asymptotically in terms of a velocity -de- 
fect layer and a wall layer. The defect layer occupies most of the boundary 
layer, and its thickness is taken equal to a boundary-layer thickness 6. 

The velocity differs from the external-flow velocity by an amount of order 

, 2 

u , the shear stress is t = 0(u ), and the layer thickness is 6 = 0(u ). 

T T T 

The much thinner wall layer has thickness denoted by 6, and the velocity 
there is small, of order u . Coordinates measured in terms of these non- 

T 

dimensional thicknesses are defined by 


y 


Y 

6 ' 


y 



6 = 0(u ) 

T 

u T — U 
~ . w.2 el 

e e T 


( 2 . 2 ) 

( 2 . 3 ) 


where 6 « 6 , and 6 has been set equal to the ratio of the nondimensional 

x/ 2 X / 2 

local kinematic viscosity and a friction velocity u (T / T) =(t/p ) ' 

^ ■'T w e w w 

based on the density at the wall. The Reynolds number Re is based on the 
geometric reference length, undisturbed external-flow velocity, and 
kinematic viscosity; all parameters are understood to be evaluated imme- 
diately upstream of the shock wave. 
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The velocity in the undisturbed boundary layer just ahead of the 
shock wave is expressed in the defect layer in terms of y and in the wall 
layer in terms of y, as follows: 



11 

0 

(2.4) 


y = 0 ( 1 ) 

(2.5) 


The form of the profile is shown in Fig. 1 for u^ « e « 1. Equations 

(2. 4) and (2. 5) are [l3] , respectively, the "law of the wake" and the "law 

of the wall," written here for a compressible boundary layer, and are 

taken to be asymptotic representations valid as u^ -*■ 0, with y and y held 

fixed respectively. Throughout most of the analysis also e 0 such that 

u^/ € -► 0. In the wall layer the Reynolds stress and the viscous stress 

Z 

are both of the same order as the wall shear stress t = 0(u ). Since 

w T 

Y = 0(6) is extremely small, the momentum equation gives t ~ As 

7= Y/6 -*■ oo, the viscous stress becomes extremely small, while t re- 
mains equal to in the limit, provided that also y = Y/ 6 0. The mix- 

ing length approximation fc p (ydU^/dy) = + . . . is introduced here for 

y « 1 and y » 1, where k is the von KS.rman constant, taken equal to 

0. 41. For a perfect gas with uniform total enthalpy, p T = and 

1 1 2 ~ 

T = — (y + 1) - ^Y - 1)U . Integration gives, for y « 1 and y » 1, 

1 

U = r sin{r"^(T /T u (k"^ In y + c)} (2.6) 

u w'eT'' ' ' 

l/2 1/2 

where c = constant and F = (y + 1) ' /(y - 1) ' . This i&. van Driest' s 

[22] result, with the added simplifying assumption of uniform total enthalpy. 
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Expansions of Eqn. (2, 6) for -»■ 1 + € and for 0 should agree. 


respectively, with expansions of the defect-layer velocity (2,4) as y -► 0 and 
of the wall“la;yer velocity (2. 5) as y oo. For 1 + € and -► 0, 

respectively, Fqn. (2.6) gives 

U ~ 1 + 6 + (u /K)(ln y - 2II ) (2. 7) 

U -u (T /T Iny + c) (2.8) 

u T w e 


where II is Coles' [13] profile parameter; c ~ 5. 0 and, for zero 

pressure gradient, II « 0. 5 or perhaps a little larger. Since y = (5/5)y, 
comparison of Eqns. (2. 6) and (2.7) gives 

K '^10(6/6) = (T^/T^)^/^U.(0 - u^(2nK"'- + c) (2.9) 

-1 -1 

where U^(e ) = F sin (F expansions (2.7) Euid (2.8) require, 

respectively, y 0 slowly and y -*• oo slowly as u^ 0; since 

u^ = 0(l/ln Re), from Eqn. (2.9), one might take, e.g., y= 0(u^ 

and y = 0(u ^) as u -► 0, where m > 0 and n > 0. A difference from the 

T T 

incompressible case arises because Eqn. (2,8) with y = (6/6)y does not 
agree with Eqn. (2,7). That is, the expansion as y oo of the wall-layer 
solution does not agree with the expansion as y -► 0 of the defect-layer solu- 
tion, Thus these solutions have no common domain of validity and cannot 
be matched. This type of problem has been discussed in detail by Lager- 
strom and Casten [23] , with a model example related to flow at low Reynolds 
number. In the present case, the density has different values for y = 0(1) 
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and for y = 0(1), and the difficulty is resolved by use of the solution (2. 6) 
for 6 « Y « 6 • this feature was also noted by Adamson and Feo [3] and 
by Melnik and Grossman [4] . 

The defect layer, where y = 0(1), has nearly constatnt density and is 
described in a first approximation by incompressible -flow equations. The 
domain of validity of Eqn, (2. 6) can be made to include y = 0(1) if 
K ^(In y - 211) is replaced by where the same function as 

for incompressible flow. Then 

IT = r sin{sin"^(r“^U ) + r”^(T /T)^/^u u_,(y)} (2.10) 

u e w' e T 01 

Expansion for u_^ -► 0 gives Eqn. (2.4) if y is held fixed, Eqn. (2.7) if 

y -► 0 sufficiently slowly that also u^ In y 0, and Eqn. (2.6) if 

y = (6/6)y 0 more rapidly, such that u^ In y is held fixed. The use of 

Eqn. (2, 10) was suggested by Maise and McDonald [24] , who showed that 

this assumed profile permits good correlation with experiniental data for 

adiabatic flat-plate boundary layers. Their interpretation of Eqn. (2.10) 

notes that a transformed velocity F sin ^(F U^) is predicted to have the 

1 /2 

incompressible form (e ) + ^t^OI^^^ everywhere outside the 

wall layer. 

A second relation between 6 and u^ for 8P/0X = 0 can be found with 
the help of the von Karman integral c£ the momentum equation, following a 
derivation similar to that for incompressible flow given, e.g. , by 
Cebeci and Smith [25] . The result is, to second order in u^/U^, 
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( 2 . 11 ) 


m. 


a U T , 

. T ,2_e,_wl/2 

^ " U u. 't ' 


oo 


m. 


(3 


■ X “5 


u. 


U 


The positive constant m is- defined by 


oo 


0 


( 2 . 12 ) 


and occurs in another context in the following section. For analytical 
purposes, the function represented in Coles' [13] forrh 


UofCy) = K ^Iny-TlK ^(1 + cos Try) (2. 13) 

for 0 < y < 1, with ^ Y ^ X* 

In the derivation which follows, the boundary-layer thickness is taken 
as one of two important characteristic lengths. The other length is the dis- 
tance from the wall to the sonic line in the undisturbed boundary layer, de- 
noted in nondimen sional form by 6 Substituting Eqn (2. 13) in Eqn. (2. 10), 
setting y = 6^/6, and expanding for 6 ^, / 5 -► 0 gives 


a k‘^1d(5/5^) = (T /T - U.(0)] - 2H k (2.14) 

As € -*■0, ln(6/6^, )~/c u^^ e [l-(\ -l)e/4+...] - 2ll; thus 6^/6 -*■ 0 

if “*■ -An alternate form of the velocity profile (2. 6) in terms of a 

* / 

coordinate y =Y/6^, is 

U = r sin{sin“^(r"^) + r"^(T /T )^/^K ‘^u In y"^} (2.15) 

u w'e T 

for y -► 0 cind y -*■ oo. 
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3. Interaction Along a Plane Wall 


As u ^0, the orders of magnitude of the mean pressure gradient and 

T 

fluid acceleration near the shock wave are larger than in the undisturbed 
boundary layer. The Reynoldsrstress transport equations can be used to 
show that in most of the boundary layer the contributions to the mean forces 
resulting from changes in the turbulent stresses are sufficiently small, in 
comparison with the pressure and inertia terms, that they may be neglected 
as u^ 0, not only in a first approximation but also in the calculation of 
some higher-order terms. Correct asymptotic representations of the mean 
velocity and pressure perturbations can therefore be derived using inviscid- 
flow equations. Also, as noted at the end of this section, displacement 
effects resulting from flow changes very close to the wall are extremely 
small, and so the largest terms in the solution for V should approach 
zero as the distance from the wall decreases. 

In the equations which follow, all laminar and turbulent stresses are 
neglected, as are the entropy changes across the shockwave; order-of- 
magnitude estimates given at the end of this section show that the neglected 
terms are in fact of higher order than any of the terms retained. The equa- 
tions describing the fluid motion can then be written in the following form: 
a div q = q • V (3.1) 

^ (y +1) “ - l)q^ (3.2) 

pq»Vq = --Y^VP (3.3) 
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Here q, q, and a = (P/p) are, respectively, the velocity vector, the 
magnitude of the velocity, auid the sound speed, all nondimens ional with the 
critical sound speed in the external flow just ahead of the shock wave. The 
gradient and divergence operators imply differentiation with respect to the 
nondirriensional variables X and Y. Crocco ' s theorem, simplified by the 
assumption of uniform total enthalpy, is 


^2Xq='y TVs (3. 4) 

— ^ —V 

where ^2 = curl q, and the specific entropy s has been made nondimension- 

al with the gas constant R. Since the upstream value of V contributes terms 
of higher order than those to be retained here, the shock-polar equation 
becomes 


— 


U U , - 1 
u d 


(3.5) 


2U /(y +1) - (U U - 1) 
u u d 


where the subscripts u and d here denote, respectively, values immediately 
upstream and downstream of the shock wave. Since the jump in the velocity 
vector across a shock wave is in a direction normal to the shock, the shock-., 
wave slope is 


dX V, 

s_ _ d 

dY ” U 

u d 

where the shock-wave location is denoted by X = X (Y). 

s 


(3.6) 
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If the nondimens ional friction velocity is small in comparison with 
the nondimensional shock- wave strength e , the sonic line in the imdistnrbed 
boundary layer is very close to the wall, as can be seen from Eqn. (2.14) 
and Fig. 1. That is, if u_^ 0 and u^/e 0, then also 5^^^ /6 -* 0. A com- 

plete description of the local pressure changes would require both an "outei*' 
solution, obtained by taking a limit of the equations with coordinates Y/6 and 
X/A held fixed, and an ’'inner" solution, obtained with Y/6^ and X/A^ 
fixed, for suitable choices of A and A^, . The shock wave can extend nearly 
to the wall, as shown in Fig, 2, and so the upstream influence described 
by the inner solution is very small; it is shown later that ^ = 0(u^ 
where 6^/5 = 0(exp(-/c e/u^)) from Eqn. (2.14). For the outer solution, 
therefore, can be taken equal to the undisturbed velocity (2,10) or (2.4). 
The inner solution describes perturbations about the undisturbed boundary- 
layer flow, while the outer solution describes perturbations about a different 
boundary-layer flow, downstream of the shock wave; the two solutions 
should match in a proper asymptotic sense. 

For Y = 0(6) the length scale A in the downstream direction is found 

from Eqn. (3.1) and the vorticity equation to be A = 0(b^6), where 
2 2 

b (e ) = 1 - M and M is the Mach number in the external flow behind a 
o o o 

normal shock wave. Coordinates x auid y are defined by 


X = 


X 

b 6 
o 


y = 


Y 

6 


where 

b = (v -J(2y +«£ +...} 

o 4 


(3.7) 


(3.8) 
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Since the shock wave is nearly normal, the shock-polar equation (3. 5) gives 
U =(l + €)^ + 0(u ). This result suggests that throughout the flow down- 

d T 

stream of the shock wave U should be represented as a constant value 
1-6 + . . . plus small perturbations of order u^. It is convenient to sepa- 
rate the rotational part, which can be calculated from Crocco' s theorem 
(3.4), and the irrotational part, which is to be found from the solution of 
Eqn. (3.1) satisfying the appropriate boundary conditions. In the limit as 
u^ 0 with X and y held fixed, the velocity components are then expressed 
in the form 


-1 fr) 2 

U = (1 + 6 ) + U^u^ (x, y;6 ) + u^ u^ YJe ) + . . . 


+ u^ (j)^(x, y;6) + u^({>^^(x, y;e ) + ... 


(3.9) 

(3.10) 


V/b^(e ) = u^<j)^^(x, y;e ) + u^<}) ^^(x, y;e ) + . . . 

where the functions of 6 shown will be expanded below for e 0. 

The entropy s is nearly constant along a streamline, and the equation 


of state gives P = p T, since changes in p ' T' are of higher order than 

terms to be retained here. It follows that along a streamline P T ^ ^ ^ 

-^/(■Y-l) 2 

P T ' to the order required here. Substitution of T = a from Eqn. 

e u ^ 

(3. 2) then gives the pressure as 


p/p = 1 - v(u-u ) -v^(u -:i)(u-u ) + ... (3.U) 

e u u u ' ' 

Also, Crocco’ s theorem gives ~ y ^P ds/dij;, where = P ~ ** P 

and - Uy, aJ^d so fi/P ~ along a streamline. Substitution in 

the expression for allows calculation of terms in the rotational part of U: 
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(3.12) 


= (1 + Z^e + ...)uQ^(y) 

‘^ 2 ''^ ^ S '*’’** ( 3 . 13 ) 

y 


Substitution of the representations (3.9) and (3. 10) into Eqn. (3.1) leads to 
differential equations for cj) and 


d) 

^Ixx 


+ c}) 


lyy 


<t> +6 

^2xx ^2yy 


0 


(3.14) 


-u*’'* + (l+€ )(M /b f{2 +(y + <|>, )<|> 

2x 0 0 o 1 lx 


Ixx 


+ (l+€)M^ct), (u!’'^ +2<|), ) (3.15) 

o ly ly Ixy 

Expansion of the shock-wave slope (3. 6 ) gives, after integration, the shock- 

wave location x =? x (y;u , e ) as 
s T 

1 

X = u X (y;e ) + . . . = “ (l + 7 € + . . . ){ c|> (0, y;e ) - ({) (0, 1;€ )} 

S T si 26 2 1 1 

+ . . . (3. 16 ) 


where the origin of coordinates has been chosen so that x^ = 0 at y = 1. 

Thus the shock wave is located at x = 0 in a first approximation, as implied 
in Fig. 2, and the flow properties are to be studied in the quarter-plane 
X > 0, y > 0. Boundary conditions at x = 0 are found from the shock-polar 
equation (3. 5), expanded in Taylor series about x = 0; 

<})^^(0, y) = -2{ 1 + (y -l)e + ...}uQj^(y) (3.17) 

-i>2xfo,y) = -^si(y)<t'i^(o.y) +-^ U-(y - }<t>iy(o, y) 


2 


01 


(y) + 


(3.18) 


It is also required that c{) 
2 2 

approach zero as x + y 


ly(X, 0) 

“► 00 . 




= 0 and that all disturbances 
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The limiting form of <{)^ as e -*■ 0 was first given in Ref. [8] ; here a 
term proportional to e is included. The solution is expressed in terms of a 
distribution of sources along the y-axis: 

(f) (x, y) = -^{l+(y-l)e + ...} J UQ^(r 7 )ln{x^ + (y-r?)^} 

-oo 

(3. 19) 

The extended definition UQ^(-y) = gives a potential for - oo < y < oo 

which is symmetric about y = 0 and thus satisfies the boundary condition 
2 2 

there. As x + y “*■ 0, the contribution to the complex velocity is 

u (<j>, -i^i), ) = - 2u {1 +(y -1)€ + . . . }k "^(In z - 2n ) + . . . (3.20) 

T lx ly T 

where z = x + iy. The pressure at the wall found from evaluation of Eqn. 
(3. 1 1) as y 0 is 

^w'^f 2x ‘^01^’^^^’^ 

p = + (2 -y -1)6 +...}— J 2 2 ^ 

e Ox +n 

2 

where = 1 + y {2e + (2y -1)€ + . . . } is the pressure ratio across a 

normal shock wave when the upstream speed is = 1 + e. At larger dis- 
2 . 2 

tances, as x + y oo, 

m 

*^T^‘*’lx " ^‘^ly^ = u^{l+(v-l)e + (3.22) 

where m^ is defined by Eqn. (2. 12); substitution of the approximate analyti- 
cal form (2. 13) gives m^^ = 8(1 + TI)/k . That is, the integrated effect is 
that of a concentrated source having nondimen sional volume strength per 
unit length equal to {l + . . . }m^u^6 . One -fourth of this fluid appears to be 
added to the external flow in the quadrant x > 0, y > 0, Since 
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d(pU) ~ (1 - M^)dU along a streamline downstream of the shockwave, and 
2 

1 - M ~ (y +l)e , the local increase in the boundary-layer displacement 
thickness is +1)€ m^u^6 + . . as can also be found by direct calcula- 
tion, An equivalent observation was made for e = 0{u^) by Melnik and 
Grossman [5, 6] . Perturbations in turbulent stresses contribute only a 
higher- order change locally; the present result does not include the far- 
ther displacement effect which occurs on a larger length scale as a new 
equilibrium velocity profile is approached. Finally, the shock-wave shape 
found from Eqn. (3. 16) is, for y 0, 


X (y) - X (0) = ~ {l+(y y + ,..} 

S S M /C 


(3. 23) 


and, for y -► oo. 


u m 

Xs(y) -x^(0) = {1 + (y — toy 


OO 


7 / %l(n)lnr7 dT7 + . . . } 

IT 0 


(3. 24) 

Fory = Y/6 -► oo, the shock-wave displacement continues to increase, and 
should be matched with a suitable perturbed external-flow solution evaluated 
as Y 0. 

The solution for <j >2 can be found in two parts, A particular solution 
of the differential equation (3. 15) can be made to satisfy homogeneous 
boundary conditions ~ = 0 if sources are distributed 

over the entire x, y plane with the source strength chosen to be an even 
function of both x and y. The boundary condition (3.18) at x = 0 is then 
satisfied by a distribution of sources along the y-axis, with strength taken 
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to be an even function of y so that ^ 2.y zero at y = 0, as in the 

solution for tj>^. Of special interest is the total source strength found by 

2 2 

carrying, out the integrations as x + y oo, with the help of integrations 


by parts and Eqns. (3. 14) and (3. 17). The pressure, correct to order 

y 


2 2 2 - 1/2 2 2 

u_(x +y) * asx+y-*" oo, and the second-order source strength 


are found to be 



- y { u [1 + (2y -1) e + . . . 1 m + u (1 + . . . )pa + . . . } 
T ^1 t 2 


1 X 

’ 2tt 2 2 

X +y 


CO 2 ^ Z 

m, = 2(5y +9)/ u (y)dy - 2(y +1) J (f> (x, 0)dy 

2 01 ly 


(3.25) 


(3.26) 


For a constant value of y such that y » 1, P initially decreases as x in- 
creases from zero, reaches a minimum at x = y, and then increases again. 

However, there is a small error at the shock wave x = x (y) = 0(e u In y), 

s T 

-12 2 

because the largest term in Eqn. (3. 25) is 0(e u^ In y/y ), whereas the 

correct first approximation is found from the shock-polar equation as 
1 -122 2 

- “y c u m /(2-rry) , If it is desired, the accuracy of Eqn. (3.25) can 
« X JL 

2 2 2 2 

be improved near x = x^ by addition of a term - y u^y «j> 2 ^( 0 , y)/(x + y ) 

2 

with y 4^0 (0» y) approximated by its leading terms 0(ln y) and 0(1) as 

y oo; away from x = x^ the added term is smaller than the second-order 

term originally shown. 

2 2 

As X + y "^0, the perturbation velocity becomes large, and it is 
again clear that an inner solution is required. For the choice of origin 
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shown in Eqn. (3. 16) and in Fig. Z, x^(0) / 0 and so the singularity in Eqn. 

(3. 20) is displaced from its correct location through a distance 
- x^(0) = 0(u^/e ). The-domain of validity near x = 0 can be extended slightly 
by addition of a term - 2u^/c ^(1 + . . .)ln(l ~ x^(0)/z) in Eqn. (3. 20) for 
^T^'^lx *" ^*^ly^’ accomplished formally by taking a limit as x 0 

with 6 x/u^ held fixed and then constructing a composite solution. The cor- 


rection is local, and introduces only a smaller change of order e u^/ [ z ] 

when |z/x^(0)| » 1. The modification is, however, necessary for matching 

with the inner solution. A discussion of the inner solution given in Ref. [8] 

is briefly reviewed here, in a slightly modified form. For Y = 0(6^,), the 

undisturbed velocity is = 1 + O(u^), and the differential equations show 

that changes in U along a streamline are also O(u^) in a distance 
l/2 

AX = 0(u ' 6j,)- Inner variables x* and y* and disturbance velocities 

7* *T» 

u* and V* are defined by 


(k X - b 6 


o 


(V 

T 

1/2 


K T 


u* = 


(U - 1), 


{k Ty^)^/^v 


(3.27) 


(3.28) 


l/2 1/2 

where factors (k ) * have been included for convenience. Equation 
(3. 1) and the vorticity equation are then approximated by the transonic small- 
disturbance equations with prescribed vorticity: 

u’^a u' /ax" - a v"/9y°' + ' • • = o (3.29) 
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(3. 30) 


/9x -9u/9y = "l/y + 


The shock-wave relations (3,5) and (3.6) become 
2 

# 1 « 2 * * 

Vj = ~(u - u ) (u + a ) + . . . 

d 2 u d d a 


(3.31) 



(3.32) 


dy u -a, 

' ad 

where the subscripts a and d again refer to quantities immediately upstream 
and downstream of the shock wave and the shock-wave location is given by 


* 


- x^ (y*). As X* ”*• - oo, u* approaches the undisturbed form 
u* ~ In y* ; the boundary condition at the wall is v* (x*, 0) = 0; and as 
X* “► oo, y* -*• oo the solution should agree with the outer solution evalu- 
ated for X - X (0) -*■ 0, y -► 0. 
s 

Although complete solutions for 0 =*= and v* can only be obtained numer- 
ically, the asymptotic behavior is found relatively easily upstream as 
X* - oo and downstream as x* -»• oo, y^*^ oo. As x=5= -*• - oo, the solution 
has the form 


1 ' 

u* ~ In y* + e^ £(y’^), ~k" f ' (y" ) (3.33) 

2 

where f" - (In y*)k f = 0 subject to the conditions that f ' (0) = 0 and that 

incoming disturbances be absent as y* co; the latter implies 
l/2 

f ‘ ~ - (In y^>= ) ' k f as y* -*■ oo. Numerical integration gives k = 0. 59. 
Downstream a suitable class of intermediate limits should be studied. As 
y* oo, a shock wave is present and must approach the nearly normal 
shock wave described by the outer solution. Thus, for y* oo, since 
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~ In y , Eqns. (3. 31) and (3. 32) give u (0, y ) ~ - In y . If an inter- 

4c 

mediate variable y^ = y ) is introduced, with 1« rj(u^,€ ) « 6/6 , 

* . . 
then In y ~ In r] + In where the first term is large and constant where- 

as the second term is 0(1) and variable. In each of the differential equa- 


tions (3. 29) and (3. 30) the two largest terms remain of the same order if 

=jc * * 4 : *1/2 

X = 0(T]/ln T) ) and v = 0(/ln r] ); then x and y (In y ) are of the 


5{£ ❖ ^ ^ at 

Ay n 


same order. In the limit as x -^00 and y -*■ 00 with x 
»}! ❖ 

fixed. In y In x and so also x /(] 


In y ) held 


4: !jC 

Ay 


In y ) ~ X 


In X ) . For the 


derivation of higher-order terms, not to be shown here, it is convenient to 

* Sje 

make this replacement. In this limit, then, the largest terms in u + In y 

❖ _ 1/2 * , ^ ^ J * 

and (In x ) v can be written as functions of x /(y /In x ). The solu- 
tions are easily obtained and the results for U and V finally can be rewritten 


1/2 -1 * l/?-l =ie-l 

U~1+(kT^/ ) u^lny - (/c ) u^ln{x (In x ) + y } 


(3. 34) 


V ~ {(y+1)(/c xy^)'^Q^ln x*}^/^2(k tan"^{y*(ln x* /x* } (3. 35) 

Factors (1 - ~ {(y+1)(kt/^) In x where M is the local 

Mach number, appear in the locations expected for solutions of the Prandtl- 
Glauert equation. The flow is represented by superposition of a known rota- 


tional flow and an initially unknown irrotational flow, described in terms of 

2 -1/2 

perturbation velocities U - I and (1 - M ) 'V which are linear in 

{kT^^^) ^u and are functions of variables (1-M^) - b 6x (0)1/6 . 

e T ‘ o s ■*' * 

r 2 2 

and Y/6 ^ . For a limit such that [x - x^(0)] + y 0 sufficiently slowly. 
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2 i/2 

with (1 - M ) ' y/[x-x (0)] held fixed, the largest terms obtained if Eqns. 

(3. 34) and (3. 35) are rewritten in the outer variables x and y are identical 
to the largest terms found from Eqns. (3. 9) and (3, 10), with the help of 
Eqn. (3. 12) and a modified Eqn. (3. 20) in which z is replaced by x - x (0) 

+ iy. Introduction of the inner solution thus removes the logarithmic singu- 
larity which appears in the outer solutions for the velocity and the pressure 
as X, y 0. 

In the derivation of these results, terms 9 (p U* U' )/9X, 9 (p U‘ V )/9 Y, 
etc., were omitted from the momentum equation, and therefore a correspond- 
ing set of terms was omitted in Eqn. (3.1) and in the calculation of the 
changes in vorticity. Expressions for these quantities, and therefore also 
order-of -magnitude error estimates for the solutions given above, can be 
obtained from the Reynolds-stress transport equations [25] . The equation 
for 9 (pU* U' )/9X contains, in particular, terms proportional to 

p U' U' 9 U/9 X. Ahead of the shock wave p U' U' is expected to be of the 
2 

same order as p U' V , of order u^. Relative changes at the shock wave 

have been estimated [26] to be proportional to the shock wave strength and 

2 

are therefore small. Thus, p U' U' = O(u^) downstream of the shock also; 

since 9 U/9 X = ){e”^^^) for X = O(e^^^u^), the product is 0(u^e 
Other terms involving velocity correlations are likewise at most 
0(u^e Neglected terms in the expansion of Eqn. (3. 1) and the vorticity 

T 

equation are also of this order, and can easily be shown to be small in 
comparison with any of the terms retained. Similarly, the derivative of 
the entropy along a mean streamline contains terms proportional to 
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p U’ V 9U/8 Y, etc., and therefore is snmall enough to be neglected in the 

derivations above. At the shock wave the entropy jump for Y = 0(6) con- 

3 

tains a constant term of order e and functions of y which are of order 
2 2 

e u^, .... It can then be shown that these changes are also suf- 
ficiently small that J2/P and PT remain constant along a mecin 

streamline to the order considered here. Finally, the changes in Reynolds 
stresses become important in a sublayer where the perturbation in is 
no longer negligible in comparison with the perturbation in p UU^. For 
X = 0(e ' u ), since t = 0(u ) and = 0(e ^ ), the sublayer is de- 

T T A 

2 1/2 

fined by Y = 0(u^e ). As will be shown in detail in Part II, the relative 

change in t is 0(e), and the new term in U which contributes to a displace- 
ment effect is 0(e u ). From the continuity equation it follows that the 

T 

2 2 

corresponding term in V is 0(e Thus, as y 0, the largest term 

in the outer solution for V which satisfies a nonzero boundary condition 
2 2 

is 0( e u ), smaller than any of the terms retained above. All of the 
neglected terms arising from these effects are smaller than the terms 


retained by at least a factor of order e. 
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4. Geometric Effects and Comparison with Experiment 

The theory of the preceding section leads to a limiting form for the 
pressure distribution as u -^0 and u -►0, for unseparated flow. In the 
flow past an airfoil at supercritical speed, with a shock wave terminating 
a region of supersonic flow, the additional effect of surface curvature can 
also be important in changing the pressure distribution and delaying sepa- 
ration, as discussed below and in Part II. The boundary layer might remain 

attached for up to about 1. 25, depending on the profile shape; Re may 
7 8 

be about 5 x 10 or perhaps as high as 10 ; and the flow ahead of the shock 
wave experiences a favorable pressure gradient, with magnitude which de- 
pends on the airfoil shape, so that the profile parameter H is 

smaller than 0, 5 (e, g. , Ref. [30] ). For a combination of parameters which 
is favorable with regard to requirements of the present theory, with = 

g 

1, 26, Re = 10 , and II = 0, the relative position of the sonic line is given 
by 6^ /6 =0. 10. This value would increase as or Re decreases or as 
n increases, as seen from Eqn. (2.14). Experimental results, however, 
are not yet available with detailed local pressure measurements for values 
of the parameters which correspond to such airfoil flows and which meet 
the requirements of the theory. For all available data, either the flow is 
separated or the values of the parameters are such that the sonic line is 
not close to the wall. Nonetheless, a comparison with, data from Refs, [l] 
and [2] has been carried out, and the agreement seems favorable provided 
that corrections for geometric effects are included. 
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A wall having convex longitudinal curvature is described locally by 
1 2 

Y ~ “ KX , where K « 1 if the radius of curvature is large in comparison 

with- the reference length used in the definitions of X and Y, A local solu- 
tion for the inviscid external flow near the foot of a normal shock "wave shows 
,a discontinuity in streamline curvature [27, 28] . Ahead of the shock wave 
> 0 to provide the required acceleration toward the wall; if the flow is 
irrotational, it follows that 0. The shock-wave relations give > 0 

and < 0 downstream; therefore also > 0, whereas the tangency con- 
dition at the wall requires V < 0 as Y -»■ 0, The term in the complex 
velocity which satisfies the required conditions as X, Y -► 0 has the derivative 

b - iV^f^ ~ - (4A)K In Z + iK + 0(K) (4. 1) 

O X X 

X -1 2 

for 0 < arg Z < — ir, where Z = b^ X + iY and, as before, b^(e) ~ (y+l)e. 

The largest omitted term is of order K and is real; the value depends on 
the flow description for Z = 0(1), and is known for symmetric two-dimen- 
sional or axisymmetric nozzle flows [29] . 

Terms and of order Ke ^^^u^ ^.nd Ku^, respectively 

when Z = 0(u_^), are now added to the expansions of U and V given by Eqns. 
(3.9) and (3. 10). The rotational part of U is unchanged, and reformulation 
of the boundary -value problem for the perturbation potential shows that <}>^ 
is unchanged, whereas now depends on K, through nonlinear terms in 

the potential equation; that is, The new terms in 

2 

contribute a change in U which is O(Ku^), smaller than terms retained pre- 
viously provided that K = 9(1). Thus, to the order considered here, for 
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u /e -► 0, a curvature correction is simply added to the earlier results. 

The new term in the pressure, written in terms of x and y, is 

^ ^ {2x In 6 + X ln(x + y ) + T y “ 2y tan ^ ^ + Ax} (4. 2) 

IT X 

where the constant A is determined only if a solution is known for the external 
flow at larger distances. 

An early careful and comprehensive experimental study was carried 
out by Ackeret, Feldmann, and Rott [1] . In Figs. 3 and 4, predicted pres- 
sures are compared with their experimental results for = 1. 32, corre- 

5 

sponding to e = 0, 247, and Re = 9. 6 x 10 , based on distance to the shock 
wave. Eqns. (2.9), (2.11) and (2.14) are used for approximate evaluation 
of other parameters. One more experimental value is needed; 6^ is chosen 
since it is easily read from the measured velocity profile and since only 
In 6 enters the equations, so that an error has small effects on other quan- 
tities, For 6^ = 0.0055, the calculations give u = 0.051, 6 = 0.021, and 

T T 

n = 0.28. This value of II seems plausible (e.g. , Ref. [30] ) because of 
the observed small favorable pressure gradient ahead of the shockwave. An 
adverse gradient of about the same magnitude is evident downstream, and is 
estimated here by P^^8P/9X ss 0. 12, where P^ is the upstream stagnation 
pressure. A corresponding term is added to Eqn, (4. 2) and the term propor- 
tional to K6 Ax is neglected. The local curvature of the plate can be inferred 
from measured pressures immediately behind the shock wave. It is estima- 
ted that P ^8 P/9 Y ss 0. 15; since P„ ~ -y V , it follows that K « 0.2. With 

t JL .A. 

the kind assistance of Prof. Z. Plaskowski of the Institiit fur Aerodynamik, 
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ETH Zurich, the author was able to measure ordinates of the plate actu- 
ally used in the experiments; values in an appropriate neighborhood confirm 
the estimate K » 0.2, The origin x = 0 is chosen at the estimated position 
of the shock wave at the edge of the boundary layer, found using measured 
pressures outside the boundary layer together with Eqn. (3.24). 

The comparison in Figs. 3 and 4 shows that the curvature effect is 
comparable in importance with the boundary -layer displacement effect; 
addition of the curvature term leads to a more pronounced "shoulder" in 
the predicted wall pressure distribution. The longitudinal pressure grad- 
ient due to tunnel divergence is also seen to be important. At the plate for 
typical values of x, say 4 < x < 14, the prediction gives about 75 percent 
of the pressure drop below the value for a one -dimensional flow; outside 
the boundary layer, at Y/6 ~ 3.6, the agreement is somewhat better. It 
is found that the velocity in Eqn. (3.21) is closely approximated by (const. )/x 
for X > 2, so that Eqn. (3.25) for the pressure is adequate here, with the 
correction (4.2). Modest changes in the assumed values of the parameters 
do not have a major effect on the comparison; for example, at a given X, 
m^/x does not depend strongly on II because 6 increases if H decreases. 

The upstream exponential decay predicted by Eqn. (3. 33) is also shown in Fig. 3, 
in the form A P/P ~ u exp{k(x - x )}, with x taken equal to -14 for 

t T O O 

approximate agreement with experiment, A major difficulty with this com- 
parison is that the upstream sonic line lies at about y = = 0.26, and 

the shock wave ends at a still larger distauice from the wall, so that the 

* * 

inner region for x = 0(1), y = 0(1) is not negligibly small. At a higher 
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Reynolds number and therefore a lower u^, the shock wave would extend 
closer to the wall, and the size of the region in Fig. 3 where no prediction 
is given would be smaller. A second serious difficulty arises because the 
flow probably was separated. The authors of Ref. [l] stated that reversed 
flow would not be ascertained at any point; however, the velocity profiles 
shown seem inconclusive, since measurements were not possible very 
close to the wall. Calculations based on the theory of Part II of the present 
paper, for the parameter values given in Ref. [l] , indicate that the flow 
was in fact separated, with a very thin separation bubble having length 
equal to a few boundary-layer thicknesses. The effect of such, a bubble 
would give a more gradual pressure rise in the region of greatest disa- 
greement in Fig. 3. Finally, a slight unsteadiness in the shock wave posi- 
tion would also contribute to a decrease in the measured pressure gradient, 
A correction for flow in a circular pipe can be derived in terms of 
cylindrical coordinates x and r defined by 


+ 

X = 


X 


b (6)R R " 


^ -R - ^-R^ 


(4. 3) 


where R is the ratio of the local pipe radius •‘"o the reference length, and Y 

4 * 

is measured inward from the wall, so that r = 0 at the axis. Solutions 

+ + 

are to be found for € -► 0, u^/ € -> 0 with x , r fixed. The wall shape is 

+ 2 + 
given by r = 1 + € f(X/R) with f = 0 at X/R = 0. Velocity components U , 

4 " 4 " 4 " 

V in the x , r directions can be written with the local curvature and 


boundarydayer effects shown separately; 
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u'^ = U^°>(X/R, r'*') + | ■!: + % t r^ 


+ . . . 


(4. 4) 


V^ = V^^^(X/R, r^;e) + K b^(e)V^^\x^, r^) + “ b^(€ ) (x^, r^) + 

(4. 5) 

2 

where now K = 6 f”(0) is the wall curvature at the foot of the shock wave, 
made nondimens ional with the reciprocal of the pipe radius. The terms 
and are the terms which would be present if the effects of the 

shock wave were ignored [29] . Terms proportional to K contain the local 
curvature effect, and terms proportional to u^6 /R contain the local boun- 
dary-layer displacement effect. The latter is described in terms of a ring 
source of radius r =1 located at x =0 and having volume strength per 

unit length equal to m^u^6 + . . . ; numerical solutions for - 0(1) given 

+ 

by Melnik and Grossman [6] also include this effect. For x -► oo, the fluid 
added at the source gives an increase of “ m^u^6 /R in U , shown explicitly 
in Eqn. (4. 4) . 

l/2 

The local solutions for X/R = 0(e ' ) are found in terms of a stream 
function defined by = - r^V^^\ where i=l, 2. 

The largest terms in Eqn. (3. 1), combined with the irrotationality condition, 

lead finally to 


OO — X 

, (i)/ + _ y' (i) ^ +T \ 

vb (x,r) = zL/a e rJ.(\r) 

^ n In 

n=l 


(4.6) 
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where J (X. ) = 0 for n=l, 2, 3, , so that the wall boundary condition 

In 

1) = 0 is satisfied; also ijj -► 0 as -► oo, and boundary values 


are to be specified at x =0. To the order required, the shock-polar equa- 

( 1 ) +2 

tion reduces to the Prandtl relation, and so ij; = r (1 - r )/4 at 

+ 

X =0. The condition that the ring source gives no term of order u 5 /R 

2 ^ 

+ + ( 2 ) + + 

in U at X =0 implies i[f = - /4 at x =0. Comparison with the 

( 2 ) 

wall boundary condition shows that ijj is discontinuous at the foot of the 

*1' - *i* *f* 

shock wave x = 0, r =1; the value obtained as x -*■ 0, r -*• 1 depends on 

the direction of approach. The coefficients a^^^ can be found from the solu- 
tions of Messiter and Adamson [29] or by direct calculation, and the coef- 
(2) 

ficients a are found directly; 
n 



X J (\ ) 
no n 



m. 


2X J (X ) 
no n 


(4.7) 


for n=l, 2, 3, . . . . 

For calculation of the pressure distribution and the shock-wave shape, 

( 2 ) 

it is convenient to introduce the corresponding velocity potential 4> , which 

satisfies 9 <j3^^^/9x^ = and 9 /dr^ = and which has a loga- 

rithmic singularity at x = 0, r =1. With the help of the asymptotic formi 

for J (X r^), one can show the singular part explicitly; 
o 11 
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•J- 

m oo J (X. r ) 

(2), + +. 1 y f o ^ 

' ' (x , r ) = - ^ 1 


-X X 
n 


n=l 


h J (X ) 
o n 


LzH 


n 


, ^ hi 

7T (n + — ) (r ) 


/ X + 

-7T (n + “)x 

e cos[7T (n + — )r - j] } 


- +~T7~z ^ ^ ^ xT§ ^ ^ - 4 ? } (4.8) 

Ztt (r ) ' 

where In ? = - (tt /4){x^ + i(l - r^)}, and ]f^ indicates that the real part 

-J- 

is to be taken. As x 0 and r -►1, the largest term in the complex 

( 2 ) ( 2 ) 

velocity U' * - iV' is due to a two-dimensional source of strength m^, 

in agreement with Eqn. (3.22). The change in the boundary -layer displace - 

(2) . (2) 

ment effect is then found by subtracting the source term from U - iV 

+ 

and adding the constant term which remains as x -»■ oo. If the numerically 
small contribution of the infinite series is omitted, the corresponding 
correction to the wall pressure is 


AP 
w 

P 


e 


- y u 


6 1 

T R 2 


( 


-Stt x^/4 


1 + 


1 - e 


-TT X 


TT X 


) 


(4. 9) 


+ y 1 / 

As X = X/(b^R) 0, AP^/P^ approaches a constant value - “y 

which implies an additional second -order correction to the boundary -layer 

solution found in Section 3 for X = 0(b 6). The shock-wave shape is found 

o 

-f' 

directly from the potential; in particular, as r -^1 the displacement of 
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the shock from its intersection with the axis is found by adding the pertur- 
bation potential from Section 3 to that found here, and subtracting the com- 
mon term proportional to In Y. The result is 

AX^ = - (47r€)‘^(Y+l)^^^€^/^m,u 6 ln(R/6)+0(u 6 e (4.10) 

S It t 

In Fig. 5 a comparison is made with pressures measured in a circu- 
lar pipe by Gadd [2] , for = 1.12 and e = 0. 097. The length of an equiv- 
alent flat'plate boundary layer is not a given quantity; instead, Gadd' s 
estimated value for boundary-layer thickness is used here, along with the 
estimate 6^/6 =0.45 found from the measured velocity profile. The sonic 
line is therefore still further from the wall than in the Ackeret experiment. 
Other approximate values are calculated as u^ = 0.04, 6 = 0.02, and 
n =0.1; the Reynolds number corresponding to these values is Re = 

6 X 10^. A pressure gradient due to small divergence of the test section 
is estimated downstream by 9 P/9X = 0. 06. The effect of finite pipe 
radius is seen to be about as large as the boundary-layer displacement 
effect. The upstream exponential decay is also shown, with x^ again taken 
equal to -14. Again there is a relatively large region where no prediction 
is made and where numerical solution of the transonic small-disturbance 
equations is required. Such a solution was obtained by Melnik and Gross - 
man [4] for this case and is also shown in Fig. 5. For large x the analy- 
tical and numerical solutions differ by an arnount about equal to the correc- 
tion for the change in pipe cross-section area. 
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5., Conclading Remarks 


The interaction of a turbulent boundary layer with a weak normal shock 
wave has been described here and elsewhere [3, 4 , 8] in terms of a rational 
approximation based on systematic asymptotic expansion procedures. The 
interaction is characterized by two small parameters, a nondimen sional 
friction velocity u^ and a nondimens ional shock -wave strength and limit- 
ing forms of the local solutions can be studied as u "*-0 and € 0, For the 

case ^ “*■ GO [3] , analytical solutions indicated that separation does not 
occur; solutions for u^/e held fixed [4] , with the first approximation de- 
scribed by the transonic small-disturbance equations, gave the same result. 
If, finally, 0, it appeared that analytical solutions would be possible 

and that perhaps the onset of separation could be discussed. Solutions for 
the pressure have been obtained here and will be used in Part II for the cal- 
culation of wall shear stress and a discussion of incipient separation. 

The largest terms in the pressure, of order u^, are derived quite 

easily, and a number of higher-order effects have been added. Corrections 

of order e u^ give, e. g. , a 35% change if = 1. 25. A partial solution for 

2 

terms of order u^ shows that these terms likewise are significant, typically 

6 7 

giving changes of 25% to 50% for Re = 10 or 10 . Corrections of order 
/c u^ and u^6 /R, obtained in ^.nalytical form for a wall with longitudinal cur- 
vature and for a circular pipe respectively, are found to be numerically 
important for the tests of Refs. [1] and [2] . 
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In the solutions for these higher-order terms, the dependence on the 
parameters is of coarse shown explicitly, and the relative importance of 
different effects is therefore apparent. It is not, however, possible to ob- 
tain analytical solutions in the asymptotically small inner region which 
accounts for the upstream influence. For values of the parameters cor- 
responding to actual transonic flight conditions, it is possible for this re- 
gion to be relatively small. Experimental results, however, are not avail- 
able in this parameter range; for existing data, either the flow is separated 
or the sonic line is not close to the wall. Nonetheless, some comparisons 
with such data were attempted, and the agreement seems fairly good down- 
stream from the inner region. The predicted pressures remain somewhat 
higher than the experimental values, and the correction terms calculated 
thus far are large enough to suggest that additional higher-order terms would 
be likely to give still further improvement. 

An essential feature of the asymptotic flow description in terms of u^ 
and € is the two-layer structure of the undisturbed profile, expressed by the 
law of the wake and the law of the wall. It is this property which permits 
the calculation of interaction pressures without knowledge of changes in 
shear stresses close to the wall. In other studies [31, 32] which were not 
based on use of this profile, derivation of a sublayer solution was necessary 
before the calculation of the pressure could be completed; these studies also 
introduced a linearized formulation for the main part of the boundary layer. 
In the present asymptotic description for u^/c -*-0, sublayer effects do not 


36 



appear even among the second -order terms in the pressure. Linear equa- 
tions appear naturally as a consequence of the limiting case considered, 
and the procedure for adding higher -order terms is clear. In the formula- 
tion for u /€ held fixed [4] , again the flow details near the wall do not 

T 

influence the pressure. The differential equations obtained in the limit 
are, however, nonlinear and numerical solution is required- For u^/« -► cd 
[3] , the pressure once more is found without knowledge of chctnges in the 
wall shear stress, 

A complete asymptotic description for 0< u^/e < oo is therefore now 

available, with numerical solutions obtained as u 0 and e 0 if u /e is 

T t' 

fixed and analytical solutions if oo or 0. For accurate, 

calculations in parameter ranges of practical interest, some further ex- 
tensions appear to be needed. In the present case, as 0» the neces- 

sary condition that the sonic line be close to the wall is met for a rela- 
tively narrow range of the parameters. The solutions for wall pressure would 
be more useful if a simple curve fit was introduced for the inner region, say 
by means of a straight line tangent to the source solution downstream and to 
the exponential solution upstream. The choice ~ exponential 

term was made for agreement with experiment in Fig. 3; the results shown 
in Fig. 5 suggest that the magnitude is too large and that perhaps a more 
suitable tentative value would be x^ = -10. The present solutions also suggest 
that terms of higher order than those retained in Ref. [4] are likely to be 
important for u^/e = 0(1). In this case the curvature correction would no 
longer have a simple form, but would have to be incorporated in the numerical 
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solation through the use of modified boundary conditions. Moreover, it ap- 

2 

pears that certain terms of order u^, and possibly still other higher -order 
corrections, will also be essential for numerical accuracy, in the wall 
shear (Part II) as well as in the pressure. Finally, the local interaction 
influences the potential flow at larger distances; the manner of introducing 
corrections in the external flow deserves farther study. 
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Solutions as a -*•0 with € fixed, such that the restriction to weak 
T 

shocks is removed and the flow is considered supersonic rather than tran- 
sonic^ permit estimation of previously neglected higher-order terms in 
It is found that an outer solution obtained for X/6 and Y/5 fixed can np 
longer be matched with the inner solution derived x and y held 

fixed. Instead, it appears that a set of intermediate linaits should also be 

* * 
considered, such that u^ In y is held fixed, e. g. , such that y - 

= 0(6 6 ) for 0 < n < 1. 

For each value of u^ In y such that 0 < u^ In y < In (6 /6^), 
the velocity profile (2.15) gives a value for in the range 1 < < 1 + c. 

The flow downstream of the shock wave is then described in terms of small 


perturbations about a uniform flow with properties found from the jump 
conditions for a normal shock wave, and indicated by a subscript o, as 


follows: 



u 


(A. 1) 



1 - U 


i 

(V+1)Uq - (Y - 1) 


(A. 2) 



1 - U 
Y+1 


(A. 3) 
(A. 4) 
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The perturbation quantities can be expressed in terms of variables X/(B Y) 

and either u In y or u In x * the latter is more convenient. Variables 
T T 

cr and X are defined by 


-1 X 

cr = u K In TT — r 

T JB 0 


X = 


B Y 
o 


(A. 5) 


Here tr = 0(1) and X = 0(1), B = B (o-), and 0 < cr < cr w^here cr will be 

o o 6 o 

defined in Eqn. (A. 8) such that U (cr ) = (1 + €) If cr 0, then 

o o 

2 _l/2 -1 * 

B ~ (v+l)(K'r' ) u In X , as noted following Eqn, (3.35), and X i 
o ex 


IS pro- 


^ ^ jT * 

VI 


portional to x /(y Vln x ), the similarity variable suggested for the deriv- 
ation of Eqns. (3.34) and (3.35). The relation /c ^u In (Y/6 . ) = (t - k ^u 

T T 

implies the equivalence of u In y and u In x already noted. At the shock 


wave X = X^(cr), any quantity f has the form 

f(/c ^u In y ) ~ f(cr)-K ^u f ' (cr)ln X 

T T s 


(A. 6) 


where f may be, e. g. , U^, U^, or p In particular, at the shock 

-1 * 

wave is evaluated by replacing k u^ In y with cr in Eqn. (2.15) and then 

setting f = U in Eqn. (A. 6); therefore U (cr) is defined by 
u o 


l/U^(o-) = r sin {sin“\r"^) + r’^T^/TJ^/^cr} 


(A. 7) 


where E = ('V+1)^^^/(Y"1)^^^ s-s before. For ~ is defined 


by 


cr = K ^u ln(6/6 ) + 211 k ^u . 


(A. 8) 


and can be calculated with the help of Eqn. (2,14). It also follows that 


u' K) = -(1 + 


1x1 X. 
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For the region downstream of the shockwave, Eqns. (3.1), (3.3), 


(3. 5) and (3. 6 ) suggest expansions of the flow quantities in the form 


U 


-1 


= Uo(a-) + K » O’) + 


-1 


V = B (o-)k u V f\, cr) + . . . 

O T 1 

P = P (o') + K u. P (\ , cr) + . . . 

o T 1 


P = P q(o') + ... 

The differential equations (3.1) and (3.3) lead to the system 
U, -XV, = 

Yp U U + P, = - YP U X"^U' - X "V' 
oolj^ Ij^ ‘^oo o o 

X X 


(A. 9) 
(A. 10) 
(A. 11) 
(A. 12) 

(A. 13) 
(A. 14) 
(A. 15) 


The vorticity referred to later, is found from these equa~ 

tions as 


u P' 

n = - ( U ' + ^ ) 

kY o yP^U^ 


(A. 16) 


The solutions for and P^ must satisfy boundary conditions at the shock 
wave X = X^, and it is required that 0 as X oo. The results are 


P' . 

^ = -U-,mX -^;^^ln(X +1) 

P' , 

V = ^ (y - tan“ X ) 

1 ypo^o 2 


(A. 17) 
(A. 18) 


P, = - 7 P' In 
1 2 o .2 


(A. 19) 


X + 1 


As X 0, and Pj^ have a logarithmic form consistent with the second 
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term in the expansion (A. 6). Equation (3. 6) for the shock-wave shape gives 


ir B P' 

X ~ K ” u — Y (A. 20) 

" 2yP (l-uf) 

o o 

for 0 < cr < cr_ ; thus X = 0(u ). The slope dX /dY is 0(u^'^^/(ln or 

o s T S T I J 

0(u ) as a -*■ 0 or O' -► cr^ , respectively. Since 0 as \ -► oo, the wall 

pressure is found from Eqn. (A. 2), for 0 < cr < cr^ , as 

P 1 - 

:5— - 1 = 2v ; + “(a ) (A. 21) 

(Y+1)U^ - (Y-1) 


For X = 0(6) and Y = 0(6), variables x and y are defined by 


X 

” B (o- )6 
o 6 



(A. 22) 


where B (cr ) is the value given by Eqn. (A. 4) with U (cr^ ) = (1 + c) The 
o o o o 

differential equations give 12 ~ - 11^ - (yP aJ^d ahead of the shock 

wave 12 ~ - (u^/6 )u j^j^(y). Relative changes in S2 along a streamline are 
small except at the shock wave, and so the first approximation for 12 remains 
a function of Y, which can be calculated from the shock-wave relations, in a 
form which agrees with Eqn. (A. l6) as y 0. The flow is then described as 
a uniform flow at speed 11 ^( 0 -^) with superimposed small perturbations, the 
rotational part being shown explicitly: 


U 


1+6 




P' (<r,) 

o 6 




01 


(y) + (x, y) 

T X 


(A. 23) 


V~ Ut B^(o-j)4.y(x,y) 


(A. 24) 
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For u “► 0, it follows from Eqn. (3.1) that <}> must satisfy Laplace' s equa- 

T 

tion. The shock-wave relation U U ~ 1 gives the boundary values for <{) (0, y), 

d IX ^ 

the shock-wave location being taken as x 0 in a first approximation; also 

Z 2 

6 (x, 0) = 0 and V<f) -*■ 0 as x + y -► oo. The solution is 

‘ y 


P' (cr ) oo 2 2 1/2 

4.(x,y) = - / “oiWMx +(y-n) } dn (A.ES) 
O o O 0 - oo 


This result has the same form as the solution (3.19), which is recovered if 

2 2 

the coefficient P' /(p U ) is expanded for e -► 0. As x + y -*• oo, the com- 

o o o 

plex velocity and shock-wave shape have the same form as in Eqns. (3. 22) 
and (3. 24), respectively, but the source strength is now u 6m,P'/(2-Y P U ), 
evaluated at cr = cr^. Since P ~ - V p (cr )U ((r_)U , the wall pressure is 

6 X 0006 X ^ 


P P (cr ) + — P' (cr )a x 
W O 6 IT o 6 T 


/ 


oo 


2 ^ 2 

X +r? 


dr) 


(A. 26) 


The expansion as x oo of an inner solution derived for x = 0(1) 

1/2 -1 * 

gives P /P ~1+2 y(kT ) u Inx, which does not agree with the ex- 

® w e e T 

pension of Eqn. (A. 26) as x 0. Instead, the two results are the leading 
terms of Eqn. (A. 21) expanded for cr -► 0 and cr cr^ respectively. Thus 
one can construct inner and outer limit-process solutions, but they cannot be 
matched in the usual way, and the expansions (A. 9) through (A. 12) are not 
limit-process expansions. This asymptotic structure is closely analogous 
to, and is a direct consequence of, that for the velocity in the undisturbed 
boundary layer. A solution for P^ which is valid for 0 < cr_< cr^ can be 
written in a form analogous to Eqn. (2.10), by use of a suitable correction 
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to the argument of the function P ; 


P - P (tr^ + - u X f 
w o o ir T 


oo 

~ 2 
0 X + rj 


(A. 27) 


If u -*-0 with X held fixed, the first two terms in the expcinsion of Eqn. 

T 

(A. 27) are the terms shown in Eqn, (A. 26). If u ^0 and x 0 with 

T 

u In X held fixed, Eqn. (A. 21) is recovered. The result (A. 27) would 

T 

replace Eqn. (3. 21) if e were not small. However, 0 < e < 0. 3 if 
1 < < 1.4, and even at e = 0. 3 these two results differ only slightly; 

e. g. , P^(cTg) differs from its three4:erm representation as e -*■ 0 by less 
than three percent at = 1,4. Corrections for e = 0(1) are therefore of 
far less importance than the higher-order effects described in Sections 3 
and 4. 
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Appendix B; Intermediate Solutions 


For « e « 1, intermediate solutions can be derived for « Y « 6 
and 6^^ « X « e 6. In applications where is a much smaller 

number than in the examples of Section 4, these solutions can be used to give 
higher-order terms in the region of the steepest gradient in wall pressure. 
From a theoretical viewpoint, the solutions can be used to show details of the 
higher-order matching of terms in the outer and inner solutions. The first- 
order inner solution for the wall pressure as x‘'“^ oo is P^/P^ 1 + O(u^lnx ), 
whereas the outer solution as x “^0 gives P^/P^ P^/P^ + 0(u^ Inx), where 

P^/P^ 1 + 0(e). That is, the inner and outer solutions for P^ describe 

small perturbations from different constant values. Matching of the solutions 
is nonetheless possible provided that 0; if, however, € = 0(1), the solutions 
no longer can be matched, as explained in Appendix A. Derivation of some 
higher-order terms for e 0 helps to verify that the matching procedure can 
be carried out to higher order. 

A shock wave in the external flow can extend into the boundary layer 
almost to the sonic line, which is located close to the wall in the present case 
for \Tj. « e. Near the sonic line the beginning of the pressure rise implies 
outgoing compression waves, in a weakly rotational mean flow described by 
the velocity-defect profile, which vdll coalesce to form an outgoing shock 
wave with very small downstream inclination from the normal to the wall. 

For Y = 0(6 j^), the undisturbed profile is f 1 + O(u^), and an 'hnner'* 
solution, discussed briefly in Section 3, describes changes in U which are 
also O(u^) along a streamline. In the limit as u^ "“O, nontrivial equations 
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are obtained for describing this initial change if X = 6^) and V = 

These are the transonic small-disturbance equations, with prescribed vorticity. 
For the second approximations, to be discussed below, an additional term must 
be retained in the expansion of Eqn. (3, 1), and a correction to the vorticity is 


needed. To second order in , for Y = 0(6;}-), the entropy s is unchanged 
across the shock wave, and also pu = 1 + O(u^), so that the streamlines are 
approximately lines Y = constant. The first terms in equation (3.4) give 
U Tds/dY, where ^ = 9V/9X - 9U/9Y, and changes along a streamline 

are related by T~^ AT - U ^ AU. Since U 1 and T ~ 1, and since 

Eqn. (3.2) with T = a gives AT - (Y - 1) AU, the change in vorticity becomes 
A^^ ^ - Y^AU. The total vorticity is then 


9V 9U 
9X " 9Y 


+ Y (U - U^) 


(B.l) 


where is found from Eqn. (2.15). A somewhat longer calculation based on 
expansion of the differential equation for the vorticity shows that this change 
in vorticity contains two contributions, one from the change in density and one. 
from the torque which acts on a fluid element because the pressure gradient 
and density gradient vectors are in different directions. 

Coordinates x”’', y *' and perturbation velocity components u , v"^ are 
defined in Eqns. (3. 27) and (3.28). The inner region Y = 0(6;j-), X = 0(u^^^^ 6^) 
sketched in Fig. 2 is to be understood as vanishingly small in the limit as 
u^/e->0. To the order needed in the derivations which follow, Eqn. (3. 1) and 
the vorticity equation become 


* * 

u u * - V « = 
X y 


-(y-p c:- 


t'f2. 

"x* + 
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V - u = 
X ‘ y 


du''' du'"' 

+ y a* u (u"' - u"") — + 

dy* ^ dy'" 


(B.3) 


where refers to u ' in the undisturbed boundary layer, found from Eqn. 

• 1/2 "1 
(2. 15), and is a constant defined by a* = (kT^' ) . Changes in the slopes 

rfc “* 1 / ^ u- 

of the characteristics dx ‘ /dy = u'“ are 0 ( 1 ) in the region where x'^ and 

y* are 0 ( 1 ), and so the outgoing compression waves presumably will first 
coalesce to form a shock wave at a finite value of y . The shock-wave shape 
will be described in the new variables by x "= ^g(y )* The shock-wave relations 
(3, 5) and (3. 6 ) become 


= — (u - u ) (u + u - U O' u + • • . ) 
d2ud duT u 


(B.4) 


dx"' 

s 


dy Ujj - u^ 


(B.5) 


where the subscripts u and d again indicate values immediately upstream and 

downstream of the shock wave. In obtaining Eqn. (B.4) use has been made 

of the anticipated result that u^ - u^, since the shock wave is nearly normal. 

The value of u'"' can depend on the solution, since for y = 0(1) the pressure 

disturbance extends ahead of the shock wave. But the outgoing compression 

1 /2 

waves which originate here, within a distance X = 0(u^ 6 ^) upstream, are 

overtaken by the shock wave when Y » 6 ,^, and so the upstream flow is the 
undisturbed boundary-layer flow when Y » 655 -. Henceforth the subscript u 
will refer to conditions in the undisturbed boundary layer. 

The discussion preceding Eqn, (3.34) suggests the use of variables t 
and \ defined by 
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t = In x '*' 


X = 


X 


(Inx”')^^^ Y (Y + 1 (or* u^ 


(B.6) 

(B.7) 


This definition of X is equivalent to that given in Eqn. (A. 5) if Inx^' is small. 
Solutions are now sought for u-j. 0 and t oo , such that t “*■ 0 and X is 
held fixed. If u"^ and v'‘' are now regarded as functions of the new variables 
t and X, the differential equations become 


. 1/2 


X "t 


U ’" (r u : + (1 - ^) \) + X V * 

*2 


= -{y-~) a-vi. 


T ^1/2 


u 


* yx -L yfi 

K ^ r "t + 


) + 


(B.8) 


= - 1 + or "' u (Y u ■' - — t + • • • ) + 

T 2 


(B.9) 


The forms for the first terms in an expansion for u'*' should correspond to 
the terms required at the shock wave, found from u^ ~ In y with y ‘ = x / (Xt ' ), 

and the forms for additional terms in u^ and v* are found using equations (B. 4), 
(B.8), and (B.9). It is then assumed that 

u'^''>' {t U q (X) + (Int)Uj^(X) + Uj^(X) + (t ^lnt)u2(X) 

+ t“^U2(X) + •••} + a’^u^{t^ Uq^^ ^(X) + (t Int) Uj^^^(X) 

+ tu^^^^X) + • • •} + • • • (B. 10) 

V* (t^^^ Vj(X) + (t Int) v^(X) + t v^(X) + • • • } 

+ Q''^u {t ^ v-^ ^(X )+•♦•}+••- (B.ll) 

T i 
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It will be seen that these assumed forms are sufficient to provide the first 

few terms of a self-consistent set of solutions to the inviscid-flow equations. 

When the assumed expansions (B. 10) and (B. 11) for t oo and u t 0 

are substituted into the differential equations (B. 8) and (B. 9)» it is easily- 

seen that the largest terms u^ and constants, the values of which 

are to be determined from the shock-wave conditions. In the new variables 

the shock-wave shape can be written as k = X (t). Substitution in the shock- 

s 

polar equation (B.4), with u“ found from Eqn, (Z.15), gives 


t^^^ V. (X ) + ♦ • • = 2 (t - t u (X )+•••} {t - ” In t 
Is Us 2 

■y _ 1 2 

- In X - — — or '*' u t -f • . . + t u^ (X ) + (In t ) u , (X ) 
s 4 T O s l s 

+ u, (X )+•••+ O''’' u t^ u \x ) + . . . -or* u t^ + • • • }^ ^ ^ 
Is T 0 s T 

Equation (B, 5) gives the shock-wave shape in the form 


dx t V. (X ) i- • • » 

s _ 1 s 

dy=" ■ ‘ - + ••• 

1 / 2 — I 

Integration then gives x'*' = 0(y'^/(lnx*) ) and so \ - 0(t ), Thus the 

s s 

functions of X^ in Eqns. (B. 12) and (B. 13) can all be written as series 

expansions about = 0» with InX^ - Int. The largest term in the square 

- 1/2 

root in Eqn. (B. 12) must be 0(t ), and so the first several terms inside 

the brackets must add to zero. In particular, it follows that 


O 




(B.13) 


u 


0 


(X) = - 1 




3 + y 

4 


(B.14) 
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(B.15) 


The differential equations for u^,, and are 


u| = 0 


^1 \ 


Xu- +v- =-l 


and it is also required that ^ 0 as \ oo . As X. 0, it is found from the 
solution for u^ that u^ InX. + const. The shock-polar equation (B. 12) then 
given u^(\^) = 1/2 and “ In X. 0 as X 0, The solutions are 




u (X) = In 

^ X^ + 1 


, v^(X) = 2tan ^ (B.16) 


^ s'; — 1 / 2 

The largest terms in u + Iny and in (Inx''") v'*', i. e. , 

^ sfc — 1 / 2 

Vj, respectively, are functions of the variable X = (x ‘'/y'^) (Inx ) and so 
have a self-similar form. Furthermore, if M is the local Mach number, 

(1 - M^) (Y+l)(l-U) and so, using Iny'*' Inx*"' for X = 0(1), 


1 - (Y + 1) a*u Inx'*' 

T 

2 - 1/2 

X {1 - M ) ^ X/Y 


(B.17) 

(B.18) 


2 1/2 

That is, the Prandtl-Glauert factor (1 - M ) appears in a familiar position 

in the variable X, and also appears in the velocity component V, which can be 

written V ^ (1 - or*u v, (X). Also, the solutions for u_, u, , and u^^^ 

T 1 0 10 

» c n't 3 "I" Y s'-" 2 ^ 

are consistent with the result u'*' - Iny + — or'*' u In y which would be 

d ■'4 T ' 

found if the shock wave were exactly normal. The expansions (B. 10) and 

(B, 11) for u"" and v'*' then describe small perturbations about a parallel flow 

y- * 3 + Y * y< 2 

with U<^l-a*'u Iny + — “ — (a u lny‘) , and the largest terms satisfy 
equations which are essentially the Prandtl-Glauert equations, but with 
specified non-zero vorticity, and written in a self-similar form. These 
terms are given in Eqns. (3. 34) and (3. 35). 
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The pressure is then found from equation {3. 11), and for X 


oo 


P /P ~l + 27or*u Inx* (B.19) 

w e T \ 

where P^ 1 - "Ye + • • • . The result (B. 19) was given by Adamson and 
Messiter [8], and, with the replacement lny'^<^ Inx'*", the first-order solutions 
(3, 34) and (3. 35) also are easily shown to be equivalent to their results for the 
velocity components. In addition, the first-order shock-wave shape, for 
6*^ « Y « 6, follows from equation (3.6) as 


X (Y) ^ ^ IT (lnx*)"^(Y+ 1)^^^ (or*u lnx*)^^^Y 

S /i T 


(B.20) 


1 :{: 2 * 

An attempt at using the term — (3 + Y) a t in u to obtain a correction 
in the wall pressure P^ leads to a difficulty because the second-order term in 
the pressure formula (3. 11) contains a factor Iny , whereas a result for X “*• oo 
is required which depends only on x'*\ If the replacement Iny*" ^ Inx* + InX is 
made in this expression, it might be anticipated that a higher-order term in 
u* would permit cancellation of the term in P which contains the factor InX. 
With this assumption the improved approximation for P^ is found to be 

/w 1 + ZVo-^u Inx* + ^ Y (3 Y - i) (a*u Inx*) +... (B.21) 

we T ^ T 

The calculation of higher-order terms below shows that the cancellation does 


actually occur, and so the assumed result (B.21) is in fact obtained for P^. 

Solutions are next derived for u^> u^^ \ and u^^^. The differential 

_(!)» 

equations (B. 8) and (B. 9) give u^^ =0 and 


5] 



"2 - ''2 = I f"! - 


X = 0 


"2 - = <i + "i)K - x’ 

- XVj‘^''=i( 5 Y+l)u| 


Xu' + =i(v| -ivi) (B.22) 


(1)' (1)' 

^ =~J (3Y+1) 


The shock“polar equation (B. 12), after substitution of the solutions (B. 16), 


and with the additional terms now needed, is 
TTt^/^+. . . + a*u t^^^ (X ) + - * ♦ 

T IS 

= 2 ^ { 2t + • • • } {~; — a'"' u t In t + — r-^ a'' u t In X + ■ 
4 T 2 T s 

+ (t ^ lnt)u, (X ) + t ^ u -5 (X ) + or'''u (tlnt)u/^\x ) 

^ s c s T Is 

+ a'‘ u t u,^^ ^(X )+••• + a'^ u t (In t + 21nX ) + 

T 1 s T s' 


(B.23) 


where the functions of Xg are to be expanded about the result (B. 20) 

X ^ TT t ^ , for t “*• 00 . Since the differential equations give u/^^ = constant 
s ^ ^ 


and Uj ~ (Y + 3 ) In X as X 0, it follows from the shock-polar equation 

(B.23) that Uj^^ (X) = - (Y+3)/4and 


U2(0) = U2(0) - 7 ^ = [u;^\x) = 


(B.24) 


In addition it is required that , v^, 0 as X-'- oo . The solutions are 


u^(X) - - 3 

2(X +1) 


v-(X) = j (tan ^X - y ; 

^ ^ 1 +X^ 


n ^ 1 , 1 , . 1 ,ir . _1 _2 37T 

u,(X) =— r In In— ; + — (— - tan X) + 

X'^ + l X'^+I X^+1 ^ ^ ® 

^ 1 ,1 n -1. X2 . 1 r*” ln(X^+l) „ 

^ X'^+l X'^+I ^ ^ X^ + 1 ^ X X^+1 

Uj =--(y+3) 

u/^\x) = - “ (Y- l)ln (X^+ 1) - |(Y+3)lnX - -j (5Y+ 1) 

^ ^ X'^ + l 

v/^^X) =-“(Y + 5)(tan‘\ - f ) -7 (5Y+1)-^^ 

^ X +1 


(B.25) 


52 



The improved solution for the wall pressure is, in terms of t = Inx’^, 

11 ^ 1 

P /P <*^1h-Yq’^u {Zt-lnt + ^^t Int - t -f • • * } 
w u T Z 8 

+ Y (q'-u^)^ {-|-(3Y - 1) t^ - “(3Y - 1) tint + ;^(5Y + l)t+ . . . } + . - - (B.26) 

The solution for u"' y^') as x^' oocan be compared with the expansion as 
X 0 of the result derived for x = 0(1), Retaining only the terms which in- 
crease in size, one finds for the velocity 


U - 1 - O''*' u Inx'*' + u Inlnx''’ + or^'u In - o-* u In (1 4-X.^) 

T 2 T T T 


Y-3 ,, ^2, 2 ,, Y + 3 * 2 ,, 

+ — : — (or u ) In X - — : — (n u ) In x In In x 
4 T 4 T 




3+Y , . 5Y+ 

r— ln\ - — 


i_J^l 


Inx^ + . . . (B. 27) 


As x’’"^ 0, the largest terms in Eqn. (B.27) agree with the largest terms in 

the expansion of Eqn. (3. 9) as x oo , found with the help of Eqns. (3, 12) and 
(3.20), For another comparison, the expansion of Eqn. (A. 27) for the 
pressure as e: 0 and x 0 is equivalent to the three terms shown in Eqn. 

(B.21). 

The order of magnitude of V as \ oo can be found from consideration 
of the region for \ » 1 in which changes in t can no longer be neglected in the 
X-momentum equation. For intermediate limits such that X -*• oo slowly, the 
inertia and pressure terms are still dominant, whereas for intermediate limits 
such that y ^ CO slowly, the equation remains 0. There must then be a 

special limit for which all these terms are retained. It is found rather easily 
that this balance of terms occurs for Y/X = O(u^). Solutions in this limit are 
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obtained in Refs. [3] and [4] and in Part II of the present work, for the cases 

e < < , e = O(u^), and e » , respectively. For the present case (Part II), 

2 ❖ 

a term of order u_ Inx in U, when substituted in an expansion of Eqn. (3. 1), 

3 

leads to the result that V = 0(u^ Inx' ), smaller than any of the terms in V 
retained above. It is this argument which implies the boundary condition V 0 
as the distance from the wall decreases, to all orders of magnitude considered 
in this and subsequent sections. 

The result (B.26) does include the correction term anticipated in Eqn. 

(B. 21) and also shows a number of additional terms. In one of the smaller 
terms a factor Iny'*' has been replaced by Inx'*' in anticipation, as before, that 
the omitted part ln\ would be cancelled in the next approximation. This minor 
difficulty appears to arise because a different approximation is needed if 
y = 0(1) and therefore ln\ = 0(t). Here u"‘' is not correctly approximated 
by -Iny'*' because the gas in this region has not passed through a shock wave, 
but instead has undergone changes described by the transonic small-disturbance 
equations in a region where x^ = 0(1). However, Py is extremely small if 
y"^ = 0(1) and x'‘* » 1, and so P^ is obtained from the above solution in the 
limit as X 00 . If the expansion had been carried out instead in terms of Iny^' 
rather than Inx'’, it would evidently have been necessary to derive additional 
"inner solutions for y'” = 0(1), because in this formulation P is expanded for 
Iny’^' CO and does not approach a function of x'” as the corresponding similarity 
variable becomes large. 
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Figure 2. Asymptotic representation of "outer" flow. 
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Figure 3. Pressure at wall with longitudinal curvature: Mg = 1.322, Re = 9-6 x 10^, K ** 0.2. 
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1. Introdaction 


When a shock wave impinges upon a wall, it penetrates the boun- 
dary layer along the surface and both the shock wave and the boundary 
layer are changed from their undisturbed states. If the boundary layer 
remains unseparated, these mutually induced changes take place in a 
small interaction region. For a turbulent boundary layer, it has been 
established [1-8] that an asymptotic description of the interaction re- 
gion requires a three layer structure. In the outermost layer, com- 
prising most of the boundary layer, pressure forces are much larger 
than forces resulting from Reynolds or viscous stresses so the govern- 
ing equations are those for an inviscid flow. For the limit process 
to be considered, the solutions for this inviscid flow region are those 
given in Part I of this paper [9] , hereafter referred to as (I). Imme- 
diately adjacent to the wall is the wall layer, in which viscous and 
Reynolds stresses dominate to lowest order. Between these two layers 
is the Reynolds stress sublayer [l] (referred to as the blending layer 
in reference [Z] ) in which momentum transfer toward the wall is carried 
out by turbulent means (Reynolds stresses); the dominant terms in the 
equation of motion are the Reynolds stress, pressure gradient, and 
inertia terms. 

This paper is concerned with the analysis of the flow in the two 
inner layers, the Reynolds stress sublayer and the wall layer, the goal 
being the calculation of the shear stress at the wall in the interaction 
region. As indicated above, the limit processes considered are those used 
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in (I). Thus, if e is equal to the nondimensional difference be- 

tween the velocity and the critical sound speed in the flow external to 
the boundary layer, and u^ is the nondimensional friction velocity, we 
consider limit processes such that u^ « e « 1, In previous analyses 
for 6 « u^ (Reference [l]) and e = O(u^) (Reference [2]) it was found 
that it was not possible to formulate an asymptotic criterion for shock 
induced separation. Here, it will be shown that even for € >> u^ 
there is no apparent asymptotic separation criterion. However, exam- 
ple calculations will be used to show that the equation derived for 
the wall shear stress may be used to predict conditions for incipient 
separation with reasonable accuracy. 

It is worthwhile reiterating the fact pointed out in (I) that for an 
unseparated boundary layer the solutions in the inviscid and inner la- 
yers are uncoupled. Because the inner layers are so thin, the change 
in pressure across them is negligible to the order retained and so the 
solution for the pressure found in the inviscid layer in the limit as 
the wall is approached is indeed the wall pressure. With this pressure 
distribution known, then, solutions in the inner layers may be found, 
leading to a relation for the wall shear stress. Thus, the unseparated 
flow case is a weak interaction problem. This is not the case for a 
laminar boundary layer and occurs for the turbulent boundary layer 
because the wall layer is so thin that the upstream influence of the 
interaction causes negligible lifting of the fluid from the wall; that is, 
to the order retained the V component of the velocity is zero, in the 
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in viscid layer, as the wall is approached. This point will be discussed 
again later. 

In order to complete the formulation of the problem in the inner 
layers, it is necessary to specify a closure condition. Here, we use 
a mixing length model, including the van Driest damping factor, to 
write an eddy viscosity [lO]. Such a closure model appears to give sat- 
isfactory results as long as the flow is unseparated [ll] and has the 
virtue of simplicity; when the flow is separated, use of such a model 
gives results which have the correct trends but which do not agree well 
with experiment. 


2. Solutions in the Inner Layers 

As in (I), transonic flow over a flat plate with a turbulent 
boundary layer is considered, with a normal shock wave intersecting 
the boundary layer; an adiabatic wall is assumed as are conditions 
such that the total enthalpy may be taken to be uniform throughout the 
flowfield. Nondimen sional Cartesian coordinates X and Y are measured 
parallel and perpendicular to the wall respectively, with the origin at 
the point where the normal shock wave intersects the boundary layer. 
Lengths are made dimensionless with respect to the distance from the 
leading edge of the flat plate to the shock impingement point, L, and 
Cartesian velocity components U and V with respect to the critical 
sound speed in the flow upstream of the shock wave and external to the 
boundary layer (hereafter referred to as the external flow), a . The 
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over bars indicate dimensional quantities. The mecin temperature, T, 
density, p, pressure, P, and viscosity coefficient, p, are referred 

- %.l. -ju 

— 

to their critical values in the external flow, e. g. , T , P , etc. We 

e e 

write the Reynolds number. Re, in the usual fashion and for convenience 

•X. 

define a Reynolds number parameter. Re , as follows 


Re 


Re 


= ( 


= ( 




) 

e 


^ 

» -r* ^ 

p a L ^ 


(la) 


(lb) 


e 

The term <p' V'>/p is included in V, where the primes denote fluc- 
tuating quantities and the bracket denotes the average value. The 

A 

friction velocity, u^, is made dimensionless with respect to a^, and 
is defined in terms of the external flow density as follows 


2 

u 

T 



^wu 1 

— — = ’ U C 

p 2 e fu 

e 


(2) 


where t is the shear stress at the wall in the undisturbed flow at 
wu 

X = L, and where c is the corresponding skin friction coefficient de- 
fined as shown. Finally, the external flow velocity and Mach number 
are written in terms of a parameter, e , as 


U = 1 + € 
e 




e 


1 - (^(U^-l) 


(3a) 


(3b) 
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where for transonic flow, e « 1, As in (I)^ the problem considered 
here is one for which u « e « 1. 

T 

In both inner layers to be considered here, the characteristic 
thickness of the region is small compared to its characteristic length. 

As a result, normal Reynolds and viscous stress terms may be neglec- 
ted compared to the corresponding shear stress terms and the trans- 
verse pressure gradient is negligible, to the order retained in the anal- 
ysis. The solutions to which these inner layer solutions must match 
in a direction normal to the wall are those solutions found in (I), 
expanded in the limit as y = Y/6 0, where 6 = 6/J_i is of the or- 
der of the bo.undary layer thickness. In the limit as x = ^/A - oo, 

where A = A/L is of the order of the extent of the interaction region, 
the solutions must match with the corresponding relations in the undis- 
turbed boundary layer. It is seen, then, that the flow problem in the 
inner layers of the interaction region is formulated as a boundary la- 
yer problem with a known pressure gradient. This also helps explain 
why an additional layer (Reynolds stress sublayer) is necessary in the 
turbulent boundary layer case. That is, in either the laminar or tur- 
bulent interaction, there is an outer layer in the interaction region 
where pressure forces dominate over shear forces, and inviscid flow 
equations hold to lower orders. Obviously, then, solutions in the outer 
layer do not satisfy the no -slip condition at the wall and a new boundary 
layer must be considered at the wall. In laminar flow, a boundary 
layer is described asymptotically by a single layer and so only one 
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so called viscous sublayer is needed in the laminar interaction. How- 
ever, a turbulent boundary layer has a two layer asymptotic structure; 
as a result, two inner layers are needed to describe this boundary 
layer in the interaction region. The Reynolds stress sublayer is the 
equivalent of the velocity defect layer, as will be seen. 

With the above remarks in mind, it is possible to write a simpli- 
fied set of governing equations in which only those terms needed in 
either of the two layers considered here are retained. They are as 
follows: 


8(gUj . 9( R.Y). ^ 0 

8X 3Y 


(4 a) 


pU 


ax 


+ p V 


aY 
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ap 

ax 
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aY 
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aY „ * aY^ 

Re 


(4b) 


^ = 0 
aY 

(4 c) 

T + (^U^ = ^ 

(4d) 

P = P T = P(X) 

(4e) 

Y Re" u 

D - { 1 - exp( - 

(4f) 


where y is the ratio of specific heats, D is the damping factor, and 

2 

K =0.41 is the von KlirmsCn constant. Since terms of order u will 

T 


be retained in the solutions, it should be pointed out that terms such 
as <p*U‘>/p in the continuity equation (4a), and < p * T’ > in the 



equation of state, (4e), which are of order u^, are not included because 
perturbations from the undisturbed flow values of each of these terms 
would be of higher than second order. Since the undisturbed flow solu- 
tions 'are considered known to second order, and we are interested only 
in the perturbations from the undisturbed flow, it is not necessary to 
include the terms in question. As mentioned previously, it is assumed 
that the wall is adiabatic, and turbulent and laminar Prandtl numbers 
are unity, so that the stagnation enthalpy is constant, as in equation 

(4d). 

As shown in (^)» for the case € /u > > 1 the distance from 
the wall to the sonic line is exponentially small compared to the thick- 
ness of the boundary layer. Since the extent of the upstream influence 
of the interaction region is ordered by the thickness of the subsonic 
region, the upstream influence is confined to a region, hereafter 
referred to as the inner region, which is exponentially small in the x 
direction compared to the main part of the interaction region, hereafter 
referred to as the outer region. That is, in the x direction, the inter- 
action region actually consists of two regions, one thin compared to 
the other; in the y direction, each of these regions is subdivided into 
the three layers mentioned previously. Following the procedure em- 
ployed in (I), the solutions in the Outer region will be shown here 
in some detail. Because the upstream influence is confined to the inner 
region, the flow entering the shock wave in the outer region is simply 
the undisturbed flow at the point in question. Inner region solutions. 
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which are found using precisely the same methods employed in the outer re- 
gion, are given in reference [12] . 


Reynolds Stress Sublayer 

In the Reynolds stress sublayer, which is intermediate to the outer 
inviscid flow layer and the wall layer, inertia terms are balanced by both 
the pressure gradient and Reynolds stress terms in the equation of motion 
in the flow direction. The extent of the outer interaction region is X = 0(A), 
where, as shown in (I), 


A = b 6 6 = 0(u ) (5a, b) 

o T 

b = (y+D^'^^e - (2y+l)e/4 + ...) (5c) 

o 

If the dimensionless (referred to Li) thickness of the layer is taken to be 

2 

Y = O(^), say, then since the Reynolds stress = and from (I), 


d P/SX = 0(u^/A), the fact that the pressure gradient and Reynolds stress 

A 


terms in Eqn. (4b) must be of the same order indicates that 6 

A ^ 

Here, for convenience, we define 6, y and x as follows: 

^ = u A 

T 


0(u A). 

T 

(6a) 


Y = 



X = Ax 


(6b, c) 


The solutions to which those in the Reynolds stress sublayer must 
match as y oo are those in the outer inviscid flow layer, written in the 
limit as Y/6 = y -► 0. The equations for the U component of velocity and 
the pressure, (Eqns. (3.9), (3.12 ), (3.19 ), (3.21) and (4.2) of (I) are 
summarized here for convenience. Thus, 
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U(x, 0) = (l + € ) ^ + [ (1 + 2\ 6 +...)uQ^(y) 


+ (1 + (v -1) € + . . . ) Uj^(x) ] + . . . 

- K 6 ~ ln(C 5x) + . . . 
ir o 


P (x) 
w 

p 


0 ) 


P 

e 


1 + y (2e + (2^ - 1 )€ ^ + . . . ) 

- u 7 (1 + (27 -l)e + )u,(x) + 

T i 

+ K 6 ^ ia(C 6x)+... 

IT o 


u^(x) 


4x f 
ir J _ 


00 


(x^ -f n^) 


(7a) 


(7 b) 
(7 c) 


where P = 1 - 7 e + 0 ( c ) is the dimensionless pressure in the external 
e 

flow. The function ^Qj|^(y) describes the variation of the velocity from its 
value in the external flow in the velocity defect layer in the undisturbed 
flow. That is, if represents this velocity component, it may be ex- 
panded as 


= ""e + “t 

and ^QjL(y) variable part of the velocity distribution in the velocity 

defect layer in the corresponding incompressible boundary layer [13] . 
The form given by Coles [14] is used here 


Uoi(y) = K. \ln y - n (1 + cos Try)] 0 y < 1 


( 9 ) 


= 0 


y> 1 
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where II is Coles* profile parameter. The last terms in Eqns. 

(7a) and (7b) are due to the curvature of the wall, i. e., for a wall with 
convex curvature described locally by 

Y = - jK + . . . (10) 

where K « 1 and K -► 0 as u. -►0 and € -*■ 0 such that K/c **► 0. The 

T 

value of the constant C is found from the solution for the flow field exter- 

o 

nal to the boundary layer. 

If Eqn. (7a) is written in terms of the Reynolds stress sublayer varia- 
ble, y, the result is 

A 

U(x, 0) =(l+e) ^+u^[(l+Zye + ...)(^ln-^+{iQj^(y )+.,.) 

+ (1 + (y - 1) € + , , , ) u^ (x) ] + , , , 

- K6 “ In (C 6x) + . . , (11) 

TT O 

where 

UQjL^y^ = K- ^(In y -2II) (12) 

Thus, equation (11) is the equation to which U(x, y) must match as y oo. 

As mentioned previously, 9P/9 Y = 0 to the order retained here. This is 
easily derived from the equation of motion in the Y direction (e, g, , see 
reference [l2]). Hence, the pressure as written in equation (7b) holds 
throughout both the Reynolds stress sublayer and the thinner wall layer. 

In view of the form of equations (11) and (12), the general expansions 
for U and V are written as follows 
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U(x, y) = (1+e) ^ In -| + UQ^(y) + u^{x, y)] 


i 

T K 


A 

+ e “ ln(-| )u^^{x, y) + e u^a^^(x, y) + 


+ K6 Uj^^(x, y; 6) + . . . 


(13a) 


V(x, y) v^v^(x, y) + ... (13b) 

The corresponding expressions for the temperature, T(x, y), and density, 
p (x, y), are found by substituting equation (13a) in the energy equation, (4d) 
and substituting the resulting expression for the temperature and equation 
(7 b) in the equation of state, Eqn. (4e). If the expansions for U, V, P, and 
p and stretched vai*iables x and y are substituted into Equation (4b), the 
governing equations for u^^, found. Thus, 
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where from Eqn. (7b), Pj^, P^j^, and P^^^ are defined as follows 
A 

P^(x) = - y u^(x) 


Pll(x) = - y (2y -l)u^(x) 


(14d) 


(15a) 

(15b) 
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(15 c) 


p ( 5 J ^ ±£JC ^ 

Ic IT O 

A 

It should be noted that both u, and P, really denote tv/o terms, one of order 

Ic Ic ^ 

In 6 and one of order 1. They are written as one here for convenience. Also, 
it is found [12] that = 0(e u.^), thus confirming the result used in Part I 
that, to the order considered here, V(x, 0) = 0 in the outer inviscid flow 
layer . 

Insofar as Uj^(x, is concerned, it is seen from Equation (11) that 

Uj^ u^(x) as ^ oo. It will be shown later that the same functional dependence 

must hold as y -►0. Since P^^ = P^^(x), the solution which satisfies both 

matching conditions and the governing equation, (14a), is 

A 

P (x) 

= ^ij(x) = - = u^(x) (16) 

This result has been used in deriving equations (14b) and (14c). It iseasily 
shown [l, 12] that the solutions to equations (14b) -(14d) may be written as 
follows 








}d^ 


(17a) 


1 ^ ^ 0 


Ic 


Ic 


J ®*p{ - — 


(x-g ) 


Zk (x-^ ) 


(17b) 


}d| (17 c) 


where the B.(^) are functions to be found by matching. As y -► oo, the inte- 


gral terms in each of equations (17) go to zero exponentially and it is seen 
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that the remaining terms match with their counterparts in equation ( 11 ). 

As X -► 0, for y = constant, the solutions satisfy the shock wave jump con- 
ditions to the relevant order, as they should. As y 0 , one finds [l] the 
following asymptotic behavior for the integrals 


/ 


Bj(l) 


exp{ - 


2k(x-|) 


}de 


± 


- B . (x) ln( ) + g. (x) + . . . 
1 C.K 1 


(18 a) 


x-b B.ce) 

g.(x) = lim [ J In b] - v B.(x) (18b) 

^ n IX -5 y 1 ei 


where y = Euler' s constant = 0. 577Z1. 
e 

The solutions for U may thus be found from equation (13a), (16), and 
(17). Since, as mentioned previously, one can find the density and temper- 
ature in terms of the velocity, using the energy equation and the equations 
of state, it is seen that a complete analytical solution may be found for the 
Reynolds stress sublayer in the outer region, valid to terms of order e u^. 
It should be noted that the continuity equation could be used to find the term 
of order in V; since it is not used anywhere in this analysis, the solution 
for V is not included. Finally, it is of interest to write the solution for U 
in the limit as y -► 0 for later use in matching with the wall layer solution. 
Thus, 
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U = (1+6 ) ^ + u^[k ^In -| + 

A * 

+ € u^K ^ln(|-)[ZY - B^^(x)ln(^) + g^^(x) + . . . ] 

+ 6 u^[(y-1)Uj^(x) + Zva^^Cy) - B^^(x)ln(-^) + gj^j^(x) + . . . ] + . . . 

+ K6 [- — ln(C 6x) - (B (x)Ln 6 + B (x) ln(-^ ) + g, (x) In 6 

TT O XC iC JtC 

+ + • • • ] + . . . (19) 


Wall Bayer 

At the wall, Reynolds stresses are zero and the skin friction is, of 
course, due entirely to viscous stress. Immediately adjacent to the wall, 
then, is a layer in which, as the wall is approached, momentum transfer is 
accomplished less and less by turbulent means and more and more by 
molecular mechanisms. In this layer, Reynolds and viscous stresses 
are of the same order. The flow entering the interaction region in this 
layer has a velocity U = O(u^) and this order holds in the interaction re- 
gion as well. If the thickness of the layer is taken to be Y = 0(6), then by 

equating the orders of the Reynolds and viscous stress terms in equation 

~ * -1 ~ 

(4b), one can show that 6 = 0[(Re u^) ]. Here, in order to write 5 in 

terms of familiar quantities, we set 


= A(Re*u )"^ 

T 


A 


u T , * 

. w.l/Z Re 

IX ' T ^ ^e Re 
•^e c 


0 ( 1 ) 


(ZOa) 

(ZOb) 
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With these orders for Y and U, and since in the interaction region X = 0(A) 
and 9P/9X = 0(a /A), it is seen that, even thoagh a pressure gradient 
exists, the only terms in equation (4b) are the Reynolds and viscous stress 
terms, to the order retained. The resulting equation is easily integrated 
to give 


p D(y)(K y -X ) + ~r 

8 V ^ 


T 8 U 2 

— — = u P T 
Q T '^e W 


(21a) 


T 

W 


(x) 


T /t 
w' WU 


(21b) 


Y = 


6y 


(21c) 


where, as indicated in Eqn. (21b), the shear stress at the wall, x (x), is 

w 

made dimensionless with its value in the undisturbed boundary layer at 
X = L, so that as X - X/A - oo, x^ -► 1. Equation (21a) is essentially 
the same equation used in references [l] and [2]* the only difference lies 
in the closure conditions used. 

With the orders mentioned above for U, P, X, and Y, and for 
p = 0(1), it is easy to show [12] that V = 0(u^6/A) and to corroborate equa- 
tion (4c) to the order retained. Since U = O(u^), then from the energy equa- 

2 

tion, (4d), it is seen that T = + O(u^) and from the equation of state, 

2 

equation (4e), then, that variations in p in the Y direction are also O(u^) . 

2 

Hence, to order u^, p = p^ and as pointed out previously by Melnik and 
Grossman [2] and Adamson and Feo [l], the fact that ^ ^ q leads to the 
result that limit process expansions in the wall layer do not match with 
corresponding expansions from the Reynolds stress sublayer^ thus, this 
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difficulty arises only in compressible flow. The difficulty may be overcome 

by taking advantage of the range of validity of equation (21a). That is, in 

any intermediate limit 6 « Y « equation (21a) is still the governing 

A 

equation; it is necessary to retain additional terms only for Y = 0(6). 

Hence equation (21a) may be used to derive solutions which will match with 

A 

those found using limit process expansions in either limit, Y = 0(6) or 
Y = 0(6). Although the methods of solution used by Adamson and Feo and 
Melnik and Grossman are equivalent, the latter' s method is more straight- 
forward and will be used here. 

It is clear both from physical arguments, and from consideration of 
equation. (21a) that as Y increases such that y » 1, the viscous terms be- 
come negligible compared to the Reynolds stress terms; this is borne out 
by using the solution to be derived to compare the two terms. Also, for 
y >> 1, the damping factor D is represented by unity plus exponentially 
small terms (eqn. 5). Finally, the density may be written in terms of the 
velocity and pres sure, in general, by using equations (4e) and (4d). Since 
P = P(x) to the order retained, equation (21a) may be integrated to give, 

T T 

U(x,y) = r sin + B(x))} (22a) 

e w' e 

r = -/ ^ (22b) 

/ v-i 

where T^ = (y+l)/2 is the temperature at the wall and T^/T^ can be 
• found in terms of from the energy equation, (4d). B(x) is a function of 
integration which should be evaluated by matching equation (22a) with the 
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limit process expansion solution valid in the wall layer (y = 0(1))* that is, 
it should be found as a result of integration from the wall to the y value in 
question, using the boundary conditicn s at the wall. However, it is only 
necessary to evaluate B(x) to lowest order here, and this may be done by 
noting that if there were no shock wave, then B(x) = C, the constant from 
the undisturbed flow wall layer solution (eqn. (2.8), in (I)). Here, then, 
B(x) is written in terms of an asymptotic expansion 

B(x) = C + io(u ) B.(x) + . . . (23) 

T 1 

where w (u ) -♦ 0 as u -*■0. 

T T 

Since equations (22a) must match with equation (19), in the limit as 
y oo, ^ 0, it is seen that must have an expansion of the following 

form 

A 

2 16 
T^(x) = 1 + aj^e + a^e + . . . + u^ t^(x) + e u^ - ln(-g )Tj^^(x) 

+ € u T (x) + . . . + K6 T (x) + . . . (24) 

T li ic 

That is, as y oo, such that u^fc“^ln y = 0(1), 

A 

ulny=uln- 7 +uln— +ulny (25) 

o 

where, as shown in (I) (equation (2.9)) 


u 

T 


K 



r sin"^{U /r} -a(— +C) 

J. ® T /C 


( 26 ) 


and the expansion for shown in equation (24) follows to insure the indica- 
ted matching. Thus, if equations (Tb), (23), (2U), (25), and (26) are 
substituted into equation (22a) and the resulting equation is compared with 
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equation (19) term by term, the unknown parameters and functions in equa~ 

tions (24) and (19) may be found. The resulting solution for the shear 

stress at the wall in the outer interaction region, and the corresponding 

values for the B. and g. (calculated once the B. are known using equation (18b)) 
11 1 

in equation (19) are 


T (x) = 1 + a € + r e + . , . - u — r u (x) 

w 1 2 T 2 1 


^ a 

+ e u^x ^In (|- )(2y - a^)(y + j ) 


a 3a 

+ € u^ I u^(x)[-“ (1 - 3y - ) + Y(y + j)] 


-1 


- 2n X (2y - aj^)(v + 2 “ 


a 

1 1 


a a 

+ X ^(2y+1 - "^)(y “ “^(In X - Ygj + In 2x)| + 


2a X 

+ K 6 — In (C 6x) + . . . 

o 


(27a) 


. / Y“1 / . ^ 

1 = -^7 2 F' 


+ 2y 




= <c-"[2y+ 1 --i] 


(27 b) 
(27c) 
(27d) 


Sii = ^ - V 


(27e) 


where u^(x) is given in equation (7c). Equation (27a), then, is the solution for 
T^(x) in the interaction region, including the effects of curvature in the ex- 
ternal flow field. It has, in most respects, the same form as the equation 
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derived by Melnik and Grossman [2]^ differing mainly in the order of the various 
terms, the inclusion of specific analytical solutions at each order of approxi- 
mation, and the inclusion of the curvature terms. 

The order to make numerical calculations for a given Reynolds number. 

Re, and external flow Mach number, M = 1 + (Y+l)e /2 + . . it is necessary 

e 

to provide relations for u_^ and 6 in terms of Re and 6 . One of the required 
equations is equation (26), with equations (20) for 6; the other is given in (I), 
(eqn. (2. ll), 2.12))* This equation, with the values of the integrals as given 
by Cebeci and Smith [lO] is repeated here for completeness. 


u u _ 

c _ _L 4. /_Ln 2 

U (1+n ) "^ ^ 

e e 


-1 ^ 

(1+n )sin (— ) 


U U , 

-S(l --S)-V2 
r' ^ 2 ' 


1 1 [2 + tV^' J 

+ - ( 

^ (1+n)^ [i-{u/r)“] 


■j- (2 + 3. 1787H + |- n 


(28) 


where T^/T^ = 1 - (U^/^) from the energy equation, (4d). Finally, it is neces 
sary to write an equation for the viscosity, fJ.(T), to be used in equations (20). 
Here, |4, = t"^ was used, with calculations being performed for n = 3/4. 

Finally, it should be noted that although the solutions pre sented here are 
found to orders of approximation such that pressure gradient and inertia terms 
were not retained in the equation of motion in the wall layer, higher order 
solutions involving these terms have been investigated [l 2]. It was found that 
the first terms to involve the pressure gradient were of order 6/A in U and of 
order 6/u^A in Thus, they give very small corrections to the solutions 


presented. 
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3. Nomerical Calculations and Separation Criterion 

The variation of t with x for various values of external flow Mach num- 
w 

ber (and thus e ) and Reynolds nun±) ers representative of modern aircraft are 
shown in figure L The numerical computations were carried out using equa- 

3/4 

tions (27a) for p = T' ,\=1. 4, n= 1/2 » C = 5 and K = 0; these values 

seem to be suitable for flow over a flat plate. The curves show the general 

features found experimentally in the interaction region. That is, x (x) goes 

through a minimum, say (x ) . ; as M increases, (x ) . decreases, 

* ^ ' w mm e w mm 

while as Re increases (x ) . increases. Thus, the effect of increasing M 

w mm ^ e 

and therefore the strength of the shock wave is to decrease the value of x 

w 

everywhere in the interaction region and hence to force the flow toward sepa- 
ration; increasing Re gives the opposite effect. It should be noted that 

X 1 as X 0 because the solution shown in Eqn. (27a) is for the outer 
w 

interaction region, A solution for x^ valid in the iimer region can be written 

in terms of the solutions for the pressure perturbations in the inner inviscid 

flow region [12]; The solution is found in precisely the same manner as 

that illustrated here for the outer region and results in a solution similar to 

that given in Eqn. (27a). Finally, a composite solution for x could be writ- 

ten, using the solutions valid in the inner and outer regions. Because of the 

limit processes considered in this work (e » u^), this composite solution 

would show only a small variation in x^ for x < 0. However, since analytical 

solutions cannot be obtained for the pressure in the inner region, no solution 

for X in the inner region has been included here, 
w ^ 

It is not possible to compare the solution for x^ with experimental results 
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for a completely two dimensional iinseparated flow because none are available. 
In those cases where the flow was apparently unseparated (e. g. , references 
fl5, 16l)j T was not measured, and in more recent work, where t has been 
measured (e. g. , references [l7, 18, 19]) the flow is separated. In separated 
flow, the shock wave takes on a lambda configuration near the boundary layer 
and a relatively strong pressure gradient develops in the Y direction in the flow 
external to the boundary layer [l8]; the flow picture is quite different from 
the unseparated flow case considered here. 

Although experimental results for truly two dimensional flow are not 
available for comparison, there is one set of measurements in a tube in which 
the flow is approximately two dimensional [20]. Thus, if R is the dimensionless 
(with respect to L) radius of the tube, 6 /R = 0.04 to 0. 08j in addition, the 
changes in the core flow (external to the boundary layer) due to the rapid in- 
crease in the boundary layer displacement thickness through the interaction 
region give corrections which are asymptotically of higher order than those 
retained in Eqn. (27a). In presenting the tube data, Gadd fitted power law 
velocity profiles to the measured profiles and inferred values of 6 (dimen- 
sional boundary layer thickness) immediately upstream of the interaction. 

Using equations derived using power law profiles, he also gave Reynolds num- 
bers associated with the tunnel stagnation pressure and Mach number for 
each test. Skin friction measurements were derived from Stanton tube 
measurements. In determining the flow parameters to be used in calculating 
for comparison with each of Gadd' s experiments, it was decided to use 
Gadd' s values of 6, Re, and stagnation pressure, as being a self- 
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consistent set of data to calculate the necessary e , 6, and II for use in 

T . Thus, it is easy to show that if u = T^, 
w 



Here, a and v are the dimensional speed of sound and kinematic viscosity 
w w ^ 

respectively, evaluated at the wall temperature (atmospheric temperature, 

taken to be 59^^.) and in the case of v^, at atmospheric sea level density; 

P refers to atmospheric pressure. Using the given values of M , P, 6, and 
at e t 

T^, Eqns. (3), (26), (with Eqn. (20) for 6), (28) and (29) were used to calcu- 
late the equivalent e , u^, 6, and II . From equation (2), then, the correspond - 
iner c^ could be calculated. The ratio of the calculated C_ to the value in- 
ferred from the standard tube measurements ranged from 1.05 to 1. 29 in 
four cases reported by Gadd (Figures 25 to 28, reference [20]). For this 


reason and because of uncertainties in the calibration of the Stanton tube, it 


was decided to compare values of c^/ which is equal to as given in 
equation (27a). The results of this comparison are shown in figure 2, for the 
case = 1.15,, Re = 7 x 10^, (Figure 25, reference [20]). The remaining 
parameters are given under figure (2). The point X* / 6 =0, defined by Gadd 
as the position at which ~ 528, was found by using equation (7b). 

It is seen that the measured upstream influence is not small. That is, 
c.^/ is not small for X < - 1 say, as required for this theory, so that even 
if the solution for for x < 0 were available, it is not expected that it would 
give good agreement. In fact, using the above mentioned parameters. 
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6 , /6 = 0.5, where 5^^ is the dimensionless distance to the sonic line in the 

5^ ' •v 

undisturbed boundary layer; evidently the values for Re, 6 , and II do not form 
a good combination for comparison with the theory. On the other hand, a 
slight unsteadiness in the position of the shock wave could have contributed 
to the slow variation of the measured upstream of and in the neighbor- 

hood of the minimum. Nevertheless, the value and the position of the minimum 
of ^ £ are predicted quite accurately. Downstream of the minimum the 
comparison is fairly good; in this regard, however, it is interesting to note 
that the negative curvature seen on the calculated curve but not on this partic- 
ular experimental curve, is a feature found in other experimental results 
which could not be used here because small separation bubbles existed. 

It is of interest at this point to consider the problem of predicting condi- 
tions under which the interaction brings the flow to the point of incipient sepa- 
ration. First, it is seen from equation (27a) that there is no asymptotic condi- 
tion for incipient separation; that is, unlike the laminar case, in which 

-1/5 

e ^ = 0(Re ^ ) is the asymptotic criterion [2l]^thereis no relation between e 
and Re which holds in the limit as Re -► oo as a condition for separation. 

This is an important difference between the two flows, and it is of interest to 
investigate the reason for its occurrence. The effect of the interaction, 
through the induced adverse pressure gradient, is to slow the fluid, in the 
boundary layer, then, the stream tubes must become wider and, due to the 
constraint of the wall, the V velocity component increases at points away 
from the wall, causing the outer flow to lift away from the wall also. In the 
laminar case, the thickening of streamtubes is greatest in the viscous sublayer. 
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The resulting V component of velocity is large enough that the flow external 
to the boundary layer is affected to lowest order so that the external and 
boundary layer flows must be considered simultaneously, i. e. , a strong 
interaction results [21, 22], No matter how large Re becomes, this strong 
interaction occurs, with the thickness of the viscous sublayer and boundary 
layer decreasing as Re increases, according to their asymptotic dependence 
on Re. The sublayer momentum flux and viscous stresses decrease and the 
strength of the shock wave necessary to cause enough displacement of . the fluid 
to result in separation decreases as Re increases. In the turbulent case, even 
for e >> u , the interaction is a weak interaction to lowest order be- 

T 

cause the wall layer is so thin. Thus, until separation occurs, the outward 
displacement of the fluid in the wall layer due to the interaction is too small 
to cause any effect in the lowest order solutions in the flow external to the 
boundary layer. A strong interaction does not occur until a separation bubble 
exists. Since there is no mechanism through which variations in the wall 
layer and external flows can interact, before a separation bubble is 
formed, it appears that no asymptotic criterion exists for incipient separation. 
However, it may be that such a criterion will result from an asymptotic solu- 
tion for the separated flow problem in the limit as the size of the bubble shrinks 
to zero. 

Although the solution for t (x) does not give an asymptotic 

w 

criterion for separation, there remains the possibility that conditions 
for incipient separation can be found simply by assuming that equation (27a) is 
an accurate solution for t^(x) at values of and Re near separation. It is 
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clear from figures 1 and 2 that the solution shows the correct form 
with a minimum value^ and it is possible to calculate the correspond- 
ing values of Re and M (i. e., e) for which (t ) . =0, the condition 

° e w mm 

for incipient separation. It must be emphasized that equation (27a) is not 

being used in an asymptotic sense in such a calculation; thus in order for 

to go to zero, one or more terms in the expansion must become as large 

as the first term. Instead, we consider equation (27a) as being a good approxi- 

3 2 

mation to t (x ) in a numerical sense as long as e and u (the orders of 
w ® T 

the first terms neglected) are small compared to one. 

To illustrate the use of equation (27a) for r ^ (x) to predict conditions 
for separation, we choose the remainder of Gadd' s tubes flow experi- 
ments in which was measured [20]. That is, Gadd presented four plots 
of vs. X'/6 (Figures 25-28, reference [20]), the first of which is shown 
in figure 2. In each case he also performed oil-flow experiments, which 
indicated that in one case (used in figure 2) the flow was not separated, but 
that in the three other cases, separation did occur. Although the plotted values 
of c^ did not indicate the occurrence of separation in these three cases, it 
should be noted that the values of c^^ were inferred from measurements from 
a Stanton tube aligned facing the flow; thus, accurate reverse flow measure- 
ments could not be made. Calculations of x were made for each of 

w 

these cases, using the same method for calculating the necessary param- 
eters, as mentioned in the discussion of figure 2. The resulting values 

for (r ) . = ( c Vc ^ ) . for each case are as follows: 

w mm f fu mm 
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(1) 

M = 1.27 
e 

Re = 10^ 


(t ) . = - 0.. 080 

w min 

(2) 

M = 1.26 
e 

Re = 1. 27 X 

10^ 

(t ) . = - 0. 020 

w mm 

(3) 

M =1.34 
e 

Re = 1 . 9 3 X 

10^ 

(t ) . = - 0. 344 

w min 


Thus these calcalated results indicate that in all three cases the flow is sepa- 
rated^ in agreement with the oil -flow experiments. In case (2), the extent 

of the region where t < 0, i.e,, the extent of the separation bubble, appears 

w 

to be very small; the flow is barely separated. 

If we denote by M the Mach number of the external flow at incipient 

6 S 

separation, equation (27a) may be used, with the condition that (t 
to find M ^ as a function of Re. A typical result is shown in figure (3) for 
K = 0 and II = l/2, i.e., for conditions associated with flow along a fiat 
plate. It is seen that according to this prediction, M increases as Re 
increases. This result is in agreement with measurements made by Roshko 
and Thomke [23] for supersonic flow at high Reynolds numbers. The magni- 
tude of the increase in over a large range of Re, however, is small 

enough that this result could help explain the conclusion that there was little 
or no variation with Re, reached by Settles, Bogdonoff, and Vas [24] . 

The effects of curvature on predicted by equation (27a) can also 

be compared with experimental results. Evidently, the only data available 
are those presented in figure (37) of reference [20], reproduced here as 
figure 4. The value of the coordinate along the abscissa, t, is given as 
t = 2. 6 6 /R where R is the radius of curvature; in view of equation (10), 
then, one can write t = 2.66K. Although there is no dependence on Reynolds 
number shown, it is assumed here that the range of Reynolds numbers is 
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6 7 

10 to 10 and calculated results are given for both values. Finally, on an 
airfoil with supercritical flow, the flow is accelerating up to the shock wave; 
since II depends on the pressure gradient in the undisturbed flow upstream of 
the interaction, the value of II on an airfoil will be different for different 
curvatures. For zero pressure gradient II -0.5, whereas for highly accel- 
erating flow n is smaller and can become negative [25], Therefore, at t = 0, 
n = 0. 5, (K = 0) and at t = 0. 015, it was decided to use a value of II for mod- 
erately accelerated flows, H = 0. The values of at t = 0 can be found 

from figure (3). Those at t = 0. 015, for which K = 0. 021 at Re = 10^ and 

7 

K = 0.028 at Re = 10 , were calculated, again using equation (27a). The re- 
sults are shown in figure (4). It is seen that at the conditions associated with 
flow over a flat plate (t=0) the calculated ^^g compares very well with the 
value given by the line drawn through the experimental data. On the other 
hand, at higher curvature (t = 0.015) the calculated values are considerably 
less than those found experimentally. In reference [20] there was some dis- 
cussion of the fact that criteria for separation might have been too stringent 
in the curved surface cases so that, for example, the point through which 

the drawn line passes at M =1. 31 perhaps should have been at M =1. 29. 

6 S 0 s 

If this were the case and if negative values of TI were called for, the agreement 
at t = 0.015, would be much better. 

The present results for criteria for shock induced incipient separation 
may be compared with theoretical predictions given by Bohning and Zierep [26], 
who postulated a two layer model for the interaction region and were 
able to calculate an equation for c^. Two comparisons were made, both for 
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flat plates, at Re values of 10 and 5 x 10 ; these Re are in the range of re- 
sults presented in reference [26]. At Re = 10^, the predicted values of M 

e s 

are 1. 24 and 1. 18 and at Re = 5 x 10 they are 1. 26 and 1. 30, where the 

calculated by the present method is the first, in each case. Thus, al- 

6 6 

though the two solutions give the same 3*-^ some Re between 10 and 5 x 10 , 

Bohning and Zierep' s solution shows a much greater variation of M with Re 

6 S 

than that shown here. However, the present results appear to be in closer 
agreement with experimental measurements [23, 24] for a related problem. 

The present theory could not be compared with very recent analytical results 
given by Inger [27] , who also used a two layer model, since conditions for 
incipient separation were not presented. 

Although there appears to be no asymptotic criterion for separation in 
the limit as Re -»• oo (u^ 0) and € -♦ 0, there remains the possibility that 

there exists a criterion involving a large as Re oo. Thus, it is neces- 
sary to consider the behavior of for e = 0(1). Based on the present anal- 

ysis, it is seen that for € = 0(1), u « 1, the solution for t (x) would be 
of the following form in the outer interaction region; 


T 

W 


(x) 


w 


(c ) + u T (x;g ) + 


w. 


(30) 


where t (« ) is the value which t (x) approaches far downstream of the shock 
w, w 

d 

wave. Since the lowest order solutions for the velocities would be of the same 

form in each of the layers as for 6 « 1, it is seen that, from matching 

solutions in the limit u -► 0, one would obtain 

T 

^w , U 

r sin {( p Tp ) ' sin" ( “ )} = U (31) 

d ' e 
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Here and are the values of P and U immediately downstream of the 

shock wave in the inviscid flow external to the boundary layer respectively, 

and are thus the values which P and U approach as x -»r oo in the interaction 

region. If the jump conditions across a shock wave are used to write and 

P^/P^ in terms of and these expressions are substituted into equation (31), 

one obtains an expression for t in terms of U . Thus, 

^ w , e 

d 


_ 


-1/1 V 

1 + (^)(U^ - 1) 
2 e 


sin ) 

e 

1 - (^)(U^ - 1) 


U 

sin (— ) 

2 e 



(32) 


If U = 1 + e is substituted into equation (32) and the resulting equation ex- 
e 

2 

panded for e « 1, it is found that = 1 + a^^e + 3 . (a ^ -1)6 /2 + . . ., in 

d 

agreement with the first three terms of equation (27a). 

It is clear from equation (32) that, since 1 < U < P for 1 < M < oo 

— e — e 

(see eqn. (36)), t ^0 for any M ; instead t goes through a minimum 

- 

d d 

value of 0. 512 at = 2. 55, for y = 1.4, cuid then begins to rise with in- 
creasing Hence, there is apparently no limiting value for 

Re -► oo. Moreover, since U -*■ P as M -»• oo, it is olear that t oo . 

e e w , 

d 

Recalling the definition of T^(x), one can see that this limit means that for 
high Mach number flow the shear stress far downstream of the interaction 
must be large compared to that of the undisturbed flow. This apparent 
anomaly can be explained by considering equation (21a). The density, in the 
first term, can be written as P/T through the equation of state. Now, in the 
wall layer, the temperature differs from the constant temperature of the 
wall by only higher order terms. On the otherhand, since 9 P/9 Y = 0 
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through the wall layer and Reynolds stress sublayer, P is the pressure from 

the inviscid flow layer and so varies from to P^. As co, P^/P^ -► oo 

and so from equation (21a), t -*■00 also. In general, since the Reynolds 

shear stress, - p < U’ V' >, includes the density, this result appears to be 

independent of the specific closure condition as long as < U' V > does not 

go to zero as -*■ 00 , and is another significant departure from the larninar 

case. Experimental verification of the large values of t at high Mach num- 

^d 

bers is given in measurements by Marvin, et al. [28] . 


4. Concluding Remarks 

The use of asymptotic methods of analysis results in a relatively simple 

relation for the shear stress at the wall in the interaction region. This relation 

may be used to predict conditions for incipient separation. In order to obtain 

the proper variation of t vs. x, which includes a minimum in t , it is neces- 

sary to include terms of higher order than the first approximation; evidently 

this would be the case also if one were to calculate t for the case u = 0(c ) 

w T 

[ 2 ] . 

Although the range of parameters in available experiments does not allow 
for exhaustive testing of the theory, comparisons which could be made are en- 
couraging; more accurate results should be obtained at the high Reynolds 
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numbers associated with modern transonic aircraft. 
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Figure 1. vs x for various values of Me and Re, for flow over a flat plate. 

(Eqn. (27a) with II = 0.5, K = O). 
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Figure 2. Comparison of calculated (eqn. (2?a)) and experimental (ref. 20) values of 
vs X’/6. Experimental conditions; Me = 1.15, Re_= 7 x 10 , wall 
temperature = 15°C (59°F)» P-^e ” 13T-9kPa (20 psig) , 6 = 0.305 cm 

(0,12 in.). Corresponding calculated parameters; e = 0.120, u^ = 0.03968, 

6 = 0.0l63ii, n = 0.312, K = 0. X'/6 = 0 at P^/Pte = 0-528 from eqn. (Ta). 
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INTERACTION BETWEEN A NORMAL SHOCK WAVE 


AND A TURBULENT BOUNDARY LAYER 
AT HIGH TRANSONIC SPEEDS 


Part III - Simplified Formulas for the 
Prediction of Surface Pressures 
and Skin Friction 


A. F. Messlter and T. C. Adamson, Jr 


1 . Introduction 


In Parts I and II an asymptotic description of the interaction 
between a fully turbulent boundary layer and a weak normal shock 
wave is derived in detail for a particular limiting case. Here in 
Part III simple approximate numerical calculations of the corre- 
sponding wall pressure and skin friction distributions are described. 


2. Undisturbed Boundary Layer 

The interaction is characterized by two small parameters, a 

nondimensional friction velocity u^ and a nondimensional shock- 

wave strength e. If is the undisturbed externalrflow velocity 

immediately outside the boundary layer and ahead of the shock wave, 

made nondimensional with the critical sound speed, e is defined by 

U = 1 + e (1) 

e 

If. is the corresponding external -flow Mach number. 


(1 + £)' 


yIl 

2 


( 2 ) 


The friction velocity is likewise made nondimensional with the 
critical sound speed, and is defined in terms of the skin -friction 


coefficient c^ just ahead of the interaction by 


2 1 2 
T 2 e f 


(3) 
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For simplicity an adiabatic wail is assumed and the total enthalpy 
is taken to be uniform. Thus the ratio of the wall temperature to 
the temperature in the external flow is 


T 

w 


T 


1 



(4) 


e 

where y is the ratio of specific heats, taken to be constant. For 
numerical calculations, the viscosity coefficient is represented by 
a power law: 



The Reynolds number Re is based on external-flow quantities just 
ahead of the shock wave and on the boundary -layer length L, the 
distance from the leading edge to the shock wave. 

The undisturbed mean-velocity profile used in Parts I 
and II is given as a transformed incompressible profile by 


r sin'^r'^u ) = U.(e) + (T /T u., 

u 1 w e T 01 


(y) 


where 


(6) 


u.(e) = rsm'‘(r'^u ), r = (7) 

1 e V“ 1 

Also, y = Y/5, where Y is the coordinate normal to the wall and 6 

is the boundary -layer thickness, both made nondimen sional with 

the boundary-layer length L». For some of the derivations the 
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function u, 


(y) is represented in Coles’ forna 


01 

Uoi(y) = K. ^In y - K (1 + cos it y) (8) 

where k ~ 0.41 is the von Karman constant and II is a profile 

parameter. For very large Reynolds numbers and for zero pres- 
sure gradient, II w 0. 5 or a little larger; for a favorable pres- 
sure gradient, as appears over most of an airfoil surface, the 
value is smaller. An equation defining the nondimens ional dis- 
tance Y = 5^ from the wall to the sonic line in the undisturbed 
boundary layer is found by setting = 1 in Eqn. (6): 

u^(c'^ln(6/6^) = (TyTJ^/^[U.(e) - U.(0)] 

- u^K"^n [1 + cos(tt 6^/6)] (9) 

where U.(0) = F sin ^(F ^). If € « 1, the first term on the right- 

1 /2 

hand side becomes (T /T ) * [ U.(e) - U. (0)1 ~ e[l - (y-l)e/4 ...] ; 

e w 1 1 

if u « e , then 6 « 6 . 

In the calculation procedures shown below, the parameters 
needed are M , Re, 6 , u , IT , and 6 ^ ; as already noted, 6 and 
6^ are dimensionless. Values for these quantities might be ob- 
tained in a number of different ways. If the first three were 
given, the last three can be calculated, as shown below; instead, 
all six quantities might be available from a numerical boundary - 
layer solution; or perhaps the nondimens ional displacement 
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thickness 6 or momentum thickness 6 might be among the given 

quantities. If needed, the definitions of 6 and 6 give, for small 

a /U , 
t' e 


6_ ~iiS (2_«L. 1, 

6 K ' T ' U 

e e 


% _ 2 +3.18n +1.5n _ ^ (5iL-l)(3 -7^- h ^ 


u 

e 

( 10 ) 


§. 1+h 2+3. 18n + 1. 511 ^ 

6 ~ ,c U " ,2 T " 2 

e 2/c e U 


( 11 ) 


For illustration, the parameters 6 will be regar- 

ded as given; then £qn. (9) relates 6^ to the other quantities, and 
two additional relations can be derived. As described in Part I, 
asymptotic matching of Eqn. (6) with a wall -layer solution gives 

a K "^ln(S /6) = (T /T )^^^U.(e) - u (2E + c) (12) 

T " ! • e w 1 T 

where 6 is the nondimen sional wall-layer thickness 


~ _ w ■ w . 1/2 e 

^ “ ti 'T ^ u Re 

*^0 6 T 


(13) 


and c « 5. 0 is a constant appearing in the wall-layer solution. The 
momentum-integral equation gives 
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1+n - 

o 

K 


U 

T 

u 

e 

, ^ 1/2 , ^ 2 +3.18n + l.Sn^ ) ^ 

't ^ T ’ 4(1 +n) ( tt 2 

' 1 e e ; U 

e 

(14) 

Equations (12) and (14) can be solved simultaneously for and II . 
In an iteration, perhaps an assumed value of II (say 0. 5) should be 
substituted into both right-hand sides; Equation (12) would then be 
solved for u , and that result would be substituted into Eqn. (14) 

T 

to give a new value of II ; the new II would again be substituted into 
both right-hand sides, etc. It can be seen from these equations 
that, since u is small, the choice of definition for 6 does not have 

T 

a strong effect on the solutions for u^ and for the product (1 + II ) 6 ; 
also, the largest term in the pressure given below is proportional 
to u (1 + n )6. 

T 

3. Pres sure Distribution 

The solutions of Part I are derived for values of u^ and e 
such that u « e « 1. The coordinate X, nondimen sional with L, 

T 

is measured along the wall, with origin taken here to be at the 

intersection of the shock wave and the edge of the boundary layer. 

The shock-wave shape is X = X (Y), where ]X (Y) j is small but 

s s 


no 



nonzero because of the interaction. The solutions are written in 

* 

terms of coordinates x and x defined by 


X = 


X _ X 

b 6 ~ A 

o 


(15) 


(k T X - X (0) X - X (0) 

e s s 


(Y+1) 


1/2^ 1/2 S 


( 16 ) 


2 2 

where b = 1 - M and M is the Mach number behind a normal 
o o o 


shock wave with upstream Mach number It follows that 


u 2.1 

e 


U 


Y+1 


(17) 


2 1 

and for € « 1, b = (y+ 1) ^ {l “ ■^2-y+l)€ +...}. The value of 

o A 


X (0) is found approximately from the results of Part I as 


2u b 6 

X (0) = - [1 + (y - + . ..] (2 + 1 . 59n) (18) 

S TT fC € C 

The pressure P is nondimensional with the critical value in 

the external flow; P and P. refer to values upstream and down- 

0 X 

stream of a normal shock wave at Mach number M ; P. is the 
stagnation pressure in the undisturbed external flow; and is 
the wall pressure. The solution for x = 0(1) is expressed in 
terms of source distributions and is approximated downstream, 
as X GO, bv P = P, +(const.)/x. The first term in the constant 
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factor is proportional to the increase in boundary -layer displace- 
ment thickness caused by the shock wave, as also found by Melnik 
and Grossman [1, E] . This first term, representing the effect of 
a source at the origin, is found to approximate the effect of the 
distributed sources quite accurately for x > 2. An expression for 
the smaller second term is found by using an approximation for the 
shock -wave shape to evaluate one of the integrals derived in Part I; 
the function of II appearing in this second -order source strength is 

then approximated by a linear function for 0 < II < 0, 5. The solu- 
* 

tion for x = 0(1) is needed only in a region which is extremely 

small if u << e and would require solution of the transonic small - 
T ^ 

ou 

disturbance equations. Upstream, as x - oo, these equations 

are approximated by linearization about the undisturbed profile, 

and the solution for x - oo is found to have the form 

P /P. ~ P /P, + u exp{k (x - X ) , where k ss 0. 59 and x is 

w' t e' t T ^ o o 

as yet unknown. 

The description given in the preceding paragraph would suffice 
for the idealized case of a flat plate in a uniform flow, where the 
pressure approaches constant values upstream and downstream of 
the interaction. For flow past an airfoil the boundary-layer effect 
can be expressed instead in terms of the pressure difference from 
a potential-theory pressure which continues to vary as the distance 
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from the shock wave increases. Upstream, for X < the ex- 

ponential solution is added to a potential -flow solution !P = P (X). 

w u 

Downstream, for X > X^, the source solution is added to a 

potential -flow solution = P^(X). Close to the shockwave, for 

X^ < X < X^, P^^ is approximated by a straight line drawn tangent 

to these upstream and downstream representations. The tangency 

conditions are sufficient to determine the equation of the straignt 

line and the values for X and X,. 

u d 

The functions P (X) and P,(X) include the effects of non uni - 
u d 

form external flow associated with nonzero values of 9 P/9 X cUid 
9 P/9 Y ahead of the shock wave. If P^(X) were known very accu- 
rately, the singular behavior dP^/dX = 0(K In X), where K is the 
local surface curvature, would appear for X 0, as explained in 
Part I. However, the potential-flow pressure distribution is pre- 
sumably found numerically, with the shock wave spread over a 
few mesh points, so that the singularity is obscured. The values 
of X^ and X^ are expected to lie outside the region of rapid pres- 
sure rise corresponding to the shock wave; if this is not true, 
some simple extrapolation of P^(X) and/or P^(X) should be used 
so that the smearing of the shock wave resulting from the numeri- 
cal calculation does not influence the pressures calculated here. 

A specific location must be estimated for the shock wave, so that 
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the origin for x can be defined. 


Three representations for the wall pressure are therefore pro- 
posed, one for each of three inter va.ls, as follows: 


X< X : — 

U Jr, 


P (X) k(x*-x *) 

Ur,. O 1 

— { 1 T u_^ e } 


(19) 


X < X < X^: 
u d 




X > X^: 
d 


w 


t t 

P.(X) 


+ ax + P 


■ymu 


Pf 2'^ X 


( 20 ) 


( 21 ) 


where P^ = { (■y+l)/2} = 1,893 for y = 1.4. The constants 

X^, X^, a, and p are determined below from the requirements 

that P and dP /dX be continuous. As already noted, k « 0.59. 
w w' 

•A* 

A value x = -14 was proposed in Part I for agreement with a 
particular experimental pressure distribution. Here Eqn. (19) 
is given in a slightly modified form, and the choice made is x^ = -10, 
partly to improve agreement with a second set of data and partly to 
avoid any implication that the value is accurate to two significant 
figures. Finally, 


m = 8 { 1 + (y-D e + . . .} + 120(1 + 6n )u (22) 

K T 

where the first term is the integrated first-order source strength 
given in Part I, and some simplifying approximations have been 


114 



r 


introduced in the second term, as explained above. 

The continuity requirements for and dP^/dX then lead to 
the following system of equations to be solved for X^, X^, a and p 


P,{X.) - P (X ) ^mu X 

slji- X J Ln_-^)+x(x - X ) A 

P^ d TT P ^ d' d u' P ^ 


P^(x^) ^ ^ 

p '"' ^d% 


(23) 


P (X ) P’ {X ) 

u u , A , u u . 

t * t 


>i^exp(k(x_^ = 


ymu 




P* (X ,) - P‘ (X ) 
d d u u 
+ z: A 


(24) 


ymu 




a - 




(25) 


p = . 


vmu 




(26) 


where 


X 


A x + X (0) 
* u s 


X = 

u 


(27) 
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(28) 


X. 


d 



Equations (23) and (24) are to be solved simultaneously for and 

Xj. If u << €, then A « A, the exponential function is numeri- 
d T ^ 

cally small, and |x | « x . Also P (X ) ~ P (0) and P,(X ) 

'u d uu u dd 

~ Pj(0), the values found at X = 0 by extrapolation, and terms con~ 
d 

taining ^d^^d^ are numerically small. These estimates 

can be used to give a first guess for x^, and the equations can then 
be solved iteratively in essentially the form given. 

The constants found from Eqns. (23) - (28) are to be substi- 
tuted into Eqns. (19) - (21) to give the wall pressure as a simple 
function of distance and a somewhat more complicated but known 
function of the parameters. To indicate the level of accuracy which 
might be expected, some of the curves from Figs, 3 and 5 of Part 
I are replotted here, with x^ now set equal to -10. A comparison 

with experimental pressures from Ref. [3] is shown in Fig. 1 for 

5 

M = 1.322, Re = 9.6 • 10 , 6 = 0.021, u = 0.051, H = 0.28, 
e T 

and 6 jjj/6 = 0. 26. Effects of nonuniform external flow, resulting 
from tunnel area change and longitudinal wall curvature, have been 
estimated and included in the predicted pressure. For a circular 
pipe, a comparison with experimental results from Ref. [4] and a 
numerical solution from Ref. [2] is shown in Fig. 2, for = 1.12, 
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Re = 6 X 10^, 6 =0.0Z, a^ = 0. 04, II =0.1 and 6^/6 =0.45. Again. 

the effects of a nonuniform external flow, in this case resulting from 

a finite pipe radius and a small area change, are included. For 

these comparisons a value for 8 ^ or for 8^/8 was regarded as 

given and the value of II was calculated using the approximate 

form of Eqn. (9), with taken to be small; use of the complete 

equation would change the theoretical curves slightly. 

The representation for the wall pressure given by Eqns. (19) - 

(21) is an approximation to the asymptotic solution derived in the 

limit as u^ -*■ 0, e 0, and 0* The scaling is then clearly 

correct if u^/ e is very small and can be shown to remain correct if 

u /e = 0(1). The complete problem formulation for u = 0(e) was 
T T 

given by Melnik and Grossman [l, 2] , in terms of the nonlinear 

transonic small -disturbance equations with prescribed vorticity, 

subject to suitable boundary conditions. For a given value of II a 

one -parameter (essentially u^/e) family of numerical solutions can 

be obtained; in their calculations II =0.5. Their results likewise 

have exponential decay upstream and source -like behavior down- 

* 

stream. Values for the constants k, m, and x could in principle 
be inferred from a sequence of numerical solutions as u /e 0; 
however, the value of m for a given II would differ somewhat from 
that given above because second -order terms included here do not 
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correspond exactly to terms in the solutions of Refs. [1] and [2] 

/ * 

which become of second order as u /e 0. A value of x obtained 

o 

from numerical solutions would differ from the value x = -10 

o 

inferred from the two measured pressure distributions shown in 
Figs. 1 and 2. 

The present results are quite simple in form and give good 
accuracy for the cases shown in Figs. 1 and 2. The numerical solu- 
tions of Refs, [l] and [2] , which give a correct approximation to 
the pressure for small values of u^ and e and for e/u^ = 0(1), 
require a considerably greater amount of computation and depend 
on two parameters. It would therefore seem that the present re- 
sults might be considered adequate for a practical airfoil calculation. 

Any approximate description of the interaction which is to be 
considered reliable should be based on a correct representation of 
the undisturbed mean -velocity profile. The profile is characterized 
by two very different length scales, a boundary -layer thickness and 
a much smaller viscous length, and so is represented in terms of 
the law of the wake and the law of the wall, described by Coles [5] 
for incompressible flow and extended for compressible jQow by 
Maise and McDonald [6] . Both the present work and that of Melnik 
and Grossman are based on this two -layer structure for the undis- 
turbed profile. The derivations of Refs. [7] and [8] , on the other 
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hand, do not properly represent the undisturbed profile and intro- 
duce approximations into the equations of motion in a rather arbi- 
trary way. In spite of these supposedly simplifying assumptions, 
the derivations remain complicated and do not show the dependence 
on parameters analytically. 

The present results, therefore, (a) are expressed by simple 
functions of distance with dependence on parameters shown explicit- 
ly; (b) are based on a systematic approximation to the equations of 
motion; and (c) for the cases shown in Figs. 1 and 2, reproduce 
numerical and experimental results quite well. Furthermore, by 
use of results given in Parts I and II the pressure distribution 
away from the wall and the shear stress at the wall can also 
be calculated. Several comments should be made, however, con- 
cerning values of constants and expected accuracy for 

* 

(1) The tentative value x^ = - 10 is a rough estimate based on a 
comparison with two measured pressure distributions, and perhaps 
should be improved. 

(2) Additional uncertainties appear in both theory and experiment. 
The shock-wave position is unsteady because of turbulent fluctua- 
tions in the boundary layer, and the maximum measured dP /dX 

w' 

is therefore at least slightly reduced. If even a very thin separa- 
tion bubble is present, the beginning of the pressure rise will be 
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moved somewhat upstream. In the theory, some simplifying 

approximations have been introduced here in the expressions 

derived in Part I for X (0) and for the second term in m. More- 

s 

over, omitted terms of still higher order than retained in Part I 
might be numerically important. 

(3) The origin x = 0 has been placed at the intersection of the 
shock wave with the edge of the boundary layer; the position 
estimated from numerical potential -flow solutions alone is not 
sufficiently accurate. Perturbations in the external flow resulting 
from, the local boundary-layer displacement effect must therefore 
be calculated, so that x can be measured from the perturbed shock- 
wave position. For the simplest correction method, P^ obtained 
for the undisturbed shock wave can be used in the boundary-layer 
equations for calculation of the displacement thickness as a func- 
tion of X in the interaction region. Calculation of the potential 
flow over the new equivalent body then gives the perturbed shock- 
wave shape, and the origin for P^ can then simply be shifted 
appropriately. In this calculation, the neglect of Py is not justi- 
fied; however, the error arises primarily for < X < X^, 
where the interpolation formula for P^ has been used. In a sense, 
then, this additional approximation can be regarded instead as a 
kind of interpolation for the displacement thickness, and probably 
introduces very little additional error. 
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(4) The theory was derived for unseparated flows with 

« e « 1, but the approximate version, with the straight-line 
interpolation, is being proposed for use in a broader range. For 
Fig. 1, as noted in Part II, it is believed that a very thin separa- 
tion bubble was present in the experiment. For Fig. 2, the Mach 
number is low enough that u^/ € ~ 0.4, which is certainly not 

very small. More comparisons with experimental data are needed 
to provide a better guide to the parameter ranges for which the 
approximation might remain reasonable. 
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4» Skin Friction 


Simple equations for the shear stress in the interaction region 

may be derived also, when the approximate forms for the pressure, 

given in Eqns. (19) and (21), are used in the equation for t . That 

is, in Part II, it was shown that t could be written in terms of 

w 

the perturbation in pressure at the wall, Pj^(x) = •* y Uj^(x), in the 
major (outer) part of the interaction region. It is easily shown [9] , 
that the same holds true in the initial part of the interaction region 
where upstream influence first causes changes from the undisturbed 


surface pressure and skin friction. Thus, from Ref. [9] , it is 


shown that in this (inner) region. 


a P (x ) 

T = 1 + U “V + • • • 

w T 2 y 


* * 

and so, using Eqn. (19) for P^ (x ), for X < X^, 

* 5 }: 

u a, k(x -X ) 

, . T 1 o . 

T = 1 + — e + . . . 

w 2y 


Likewise, if Eqn. (21) is used for the pressure perturbation, 

Pj(x) = - y Uj^(x) , in equation (27a) of Part II, one can show that, 

for X > Xj 
d 


w 

4. ^ 

T , + U 

wd T X 

+ eu (-B + Clnx) 

T 

Ola) 


m 

{l - ^[l -3v- 2^ + 


4ir 

(l+(Y-l)e) 

(31b) 
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^ 3 . 

B = - /C ”^ln(|- )(2'^-a^)(Y + "I - “^) + K ^[2H (Zv-a^) (V + ” - ) 


a, a 

+ (2y+i “ "^)(y -~y ^ 


-1 ^1 ^1 
C = K (2y+1 -^)(v - -^) 


(31c) 




(31d) 


In Eqn. (31c) v = 0. 57721 is Euler' s constant and k = 0.41 is the 
e 

Karman constant. In Eqn. (31a), represents the wall shear 

stress calculated using the pressure downstream, of the shock wave. 

For flow over a flat plate, is given by Eqn, (32) of Fart II; 

2 

for e « 1, T = 1 + a, e + a, (a, -1) c /2 + . . . as in Eqn. (27a). 
wd 111 

Finally, t = cVc^ where c. is the skin friction and the subscript 
w F fu f 

u refers to conditions upstream of the interaction. 

Equations (30) and (31a) are simple relations for up- 

stream of X and downstream of X,. In the intermediate region, 
u d " 

just as was done for the surface pressure, we use a straight line 
tangent to the solution given by Eqn. (30) at X^ and to the solution 
given by Eqn. (31a) at X^. The values for X^ and X^ found here 
are clearly not the same as those found for the surface pressure 
distribution. In addition, and c^^, the skin friction upstream 

and downstream of the interaction respectively, are considered 
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to be known functions of X. That is, the present calculations may- 
be considered to be a correction to the c^ distribution found using 
typical boundary layer methods, this distribution being in error 

within the interaction region. Because c_ and c_, would include 

fu fd 

the effects of wall curvature, terms involving these effects are 

not included here. The proposed equations to be used, then, are 

♦ ^ 

a k(x -X ) 

X<X c. = c. (X)[l+u ~e ^ -t...] (32a) 

u f fu *■ T 2y -* 

^ + “c- + Pc (32b) 

X > c, = c(X) + (X )[u - + € 11 (-B + Clnx)l (32c) 

df fd fuU TX T 

where ^ 


X = 


X 

A 


A. X +X (0) 
^ s 

A 


(33) 


It may be noted that as x becomes very large, the term - B -t C In x 
does not approach zero. This term evidently indicates that relzixa- 
tion to a final value of c^ takes place on a scale large compared to 
a boundary layer thickness. Because B is a large number 

- B -t C In X goes through zero for x large (e. g. , x « 200). More- 
over, since In x does not vary rapidly with large x, the contribution 
of this term is not significant for a considerable range of values of 
X about the point at which it equals zero. Hence after the term 

- B + C In X has gone to zero, it should be dropped from the 
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expression for c^. 

The equations necessary to find X^, X^, or^, and are 

obtained by equating the relevant expressions for and d c^dX 

at X and X,. They may be written as follows: 
u d 


C, (X ) - c„(X ) 
fu u fd d 


Ax c • (X ) - Ax c ^ (^j) 
u fd u d fd d 


+ a c. (X )[— - e(B + C(1 - In x,))] 
T fu U X , d 

d 


- Q /a c • (X ) + ^ c (X 
t\ fu u A^ fuu7 


■{ < (X ) - C ■ (X )) + c (X )(--^+ — )l 

I u fu d fu u fu u Z X - j 




ic (X)(l-kx -^)-xAc' (X)| (34a) 

(fuu VI fu } 


k(x -X ) 
u o 


^ [a c ■ (X ) + k c (X 
a , fu u A . f u u 

1 * 




X - d 
d 


- u A . 

a = A c » (X^) + (X ) ^ + e u — > 

c fu ' d' fu u' i 2 T X , i 


(34c) 


Pc = - + C (1 - lA x^))} 


where 
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C J = d c^/ dX. 

An illustration of the results obtained using the present 
simplified equations is shown in figure (3) where the experimental 
results are those of Gadd [4]. The numerical values of the param- 
eters are those given in figure Z, Part II. In addition, x =1. 54, 

S 

3jC 

X, = 1.37, X = - 6. 10, and t = c / c =0.85 and x = - 10. 
d u wd fd fu o 

If the calculated results in figure (3) and figure (3 of Part II are 

compared, one can see that the effect of using the approximate 

form for P, (x) is to change the location and value of (t ) 

1 w mm 

= (c^/c^^)^^^ slightly. Comparison of calculated and experimental 

values in figure (3) shows that the upstream influence is under - 
* 

predicted for = - 10 as in the surface pressure calculations. 
However, unsteadiness in the shock wave position, with oscillations 
of relatively small amplitude, could easily account for the difference 
between calculated and experimental values. 
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X 

'Figure 2. Pressure at wall of circular pipe: = 1.12, Re « 6 x 10^; 

calculated from theory of Part I, with effect of finite pipe radius 
(6/R ft 0.055) included. 
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